BI-ORTHOGONAL SYSTEMS ON THE UNIT CIRCLE, REGULAR 
SEMI-CLASSICAL WEIGHTS AND INTEGRABLE SYSTEMS - II. 



N.S. WITTE 



Abstract. We derive the Christoffel-Geronimus-Uvarov transformations of a sys- 
tem of bi-ortliogonal polynomials and associated functions on the unit circle, 
that is to say the modification of the system corresponding to a rational modi- 
fication of the weight function. In the specialisation of the weight function to 
the regular semi-classical case with an arbitrary number of regular singularities 
{zi,... ,Zm] the bi-orthogonal system is known to be isomonodromy preserving 
with respect to deformations of the singular points. If the zeros and poles of the 
Christoffel-Geronimus-Uvarov factors coincide with the singularities then we have 
the Schlesinger transformations of this isomonodromic system. Compatibility of 
the Schlesinger transformations with the other structures of the system - the recur- 
rence relations, the spectral derivatives and deformation derivatives is explicitly 
deduced. Various forms of Hirota-Miwa equations are derived for the T-functions 
or equivalently Toeplitz determinants of the system. 



1. Motivations 

The imitary group U{N) with Haar (uniform) measure has eigenvalue proba- 
bility density function (see e.g. Ull Chapter 2]) 

(1-1) (2^.Uj''-''\'' ^r--e'''^T, e,e{-n,nl 

where T = {z e C : |z| = 1}. One of the motivations of our study is to charac- 
terise averages over Lf e U{N) of class functions w{U) (i.e. symmetric fimctions 
of the eigenvalues of U only) which have the factorization property W^^i for 
{zi, . . .,zn} e Spec(!J). Such functions w{z) can be interpreted as weights. Intro- 
ducing the Fourier components {ivi}i(=x of the weight iv{z) = Y^'il-oo due to the 
well known Heine identity [48 1 

N 

1=1 

we are equivalently studying Toeplitz determinants. 

Such averages over the unitary group are ubiquitous in many applications to 
mathematical physics, in particular the gap probabilities and characteristic poly- 
nomial averages in the circular ensembles of random matrix theory fT2l ,lll1, IIT6ll , 
the spin-spin correlations of the planar Ising model |39l,f3r|, the density matrix 
of a system of impenetrable bosons on the ring [14J and probability distributions 
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for various classes of non-intersecting lattice path problems IITSl . All of these ap- 
plications can be subsumed within a generic class of weights, termed the regular 
semi-classical class and its degenerations. The simplest non-trivial example of a 
regular semi-classical weight is 



(1.3) 



|1, 0e(-7I,7T-(/)) 
11-5, 0e(7I-c/),7T] 



where n, oo = 0Ji + ico2 are complex parameters and ci.,t = e"? are complex variables. 
It has been previously shown (161 that the average lll.2l l with weight | |1.3| | can, as 
a function of t, be characterised as a x-function for the Painleve VI system. The 
examples mentioned above then are characterised by solutions to the sixth Painleve 
equation with the parameters given in Table[T] One of our primary motivations for 
our study is to deduce recurrence relations that can be used to characterise (also in 
an algorithmic sense) this class of unitary group averages. 





Interpretation 


^ 


Generating function for the probability of exactly k 
eigenvalues of a random unitary matrix with phases 
in the sector (tt - (p, n], Dyson CUE 


k 


A:-th moment of the characteristic polynomial for ran- 
dom unitary matrices 


1/4 1/4 -i/2 


Diagonal spin-spin correlations of square lattice 
Ising model 


1/2 2 


Density matrix of impenetrable Bose gas confined 
on a circle 



Table 1. Examples of the regular semi-classical weight (|1.3|l oc- 
curring in random matrix theory and condensed matter physics 



However one notices that members of this class of weights are not necessar- 
ily positive or even real valued, iv{z) + w(z) for z e T, where the bar denotes 
the complex conjugate, and consequently the Toeplitz matrices are non-Hermitian 
Wn + vj-n- Systems of orthogonal polynomials constructed with such weights de- 
fined on the unit circle are naturally split into a bi-orthogonal pair of polynomials 
(see I I2.3I I, l|2.4|l ) . Following the pioneering works on orthogonal polynomial systems 
(OPS) on the imit circle by Szego and Geronimus (see the monographs [48], |19| 
and the contemporary state of affairs in Simon 1,44.1 , 1451 ) the study of bi-orthogonal 
polynomial systems (BOPS) was initiated by Baxter f4] motivated in part by the 
desire to analyse non-hermitian Toeplitz matrices and applications in probability 
theory. In an independent development a study of systems of orthogonal Lau- 
rent polynomials was begun by Jones and Thron i34J which has its origins in the 



3 

theory of a particular type of analytic continued fraction, the T-fractions, in mul- 
tipoint Fade approximation and the strong forms of the Stieltjes and Hamburger 
moment problems. In subsequent work [SF], f25l it was realised that these orthog- 
onal Laurent polynomial systems (OLPS) are equivalent to a pair of bi-orthogonal 
polynomials, and furthermore that these were foimd to be equivalent to Baxter's 
BOPS in [41J. In fact one can construct a system of bi-orthogonal Laurent poly- 
nomials through the defining recurrence relations, with the consequence that our 
conclusions apply equally to the bi-orthogonal Laurent polynomial systems or 
BOLPS. 

A systematic study of BOPS on the unit circle and the regular semi-classical 
class was initiated in (17) where a number of basic structures were laid out - the 
well-known recurrence relation, the spectral derivatives of the system and the 
deformation derivatives with respect to the free singular points of the weight. 
Subsequently a similar program in the context of BOLPS was undertaken in (6|. 
The present study is a continuation of the program 1 17] and here the central theme 
is the rational modification of the weight and the consequences for the system and 
its structures. 

Why do we consider rational modification? The answer to this question is 
two-fold - there are direct applications of such formulae and because a given 
regular semi-classical weight should be seen as a member of a family of weights 
whose exponents differ by integer increments or decrements and such systems are 
linked by Schlesinger transformations. The historical literature has focused on the 
rational modification of weights defined on the real line or the interval and the 
formulae expressing the new orthogonal polynomials in terms of the original ones 
were derived - see the work of Uvarov [50J, [51 1 following the pioneering study of 
Christoffel [8J. There has emerged a modern period of interest in this subject, see 
for example [53|, |46|, and [7\. 

Formulae for the polynomials orthogonal with respect to weights changed by 
polynomial multiplication and division written in terms of the original system 
were extended to those defined on the unit circle in the works [20], {2T\, [26] and 
||22J . The extension of such formulae to BOPS on the unit circle was undertaken, 
however in an incomplete way, in l42l . The most general rational modification of 
a weight iv{z) on the unit circle takes the form 

(1.4) w 
(1.5) 

Whilst this contains redundancies it manifestly exhibits interesting symmetries 
and is rendered real and positive under simple conditions on the parameters. In 
addition it is this case, or more precisely the polynomial form, that arises in studies 
on the symmetrised model of last passage percolation and non-intersecting lattice 
paths, e.g. see Formulae (2.3), (3.9) and (6.12) of [15] which were first derived in 
||3l . The general rational modification l|1.51 is also fundamental in the calculation 
of the averages of ratios of characteristic polynomials (which are an example of a 
class function) with respect to the Haar measure (here the original weight is trivial, 
w{z) = 1). In this situation we are, in effect, considering the average (|1.2|l with 
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respect to a rational weight function. A contemporary interest in such averages 
arises from analytic number theory where the average l|1.2|l of lll.St with iv{z) = lis 
conjectured to reproduce the integral of ratios of the Riemann zeta function along 
the critical line [11], [10 1, [9|, in an asymptotic sense as N ^ oo. This conjecture has 
been enlarged to apply to similar averages of certain families of L-functions, which 
now involve the calculation of averages of ratios of characteristic polynomials over 
the Haar measure for the compact classical groups Sp{N), 0"'"(2N) and 0"'"(2N + 1). 



Our penultimate result gives the bi-orthogonal polynomials 'S'ur^n their associated functions H„, S* for the modified weight 
defined by jl.Si as a {K+ K* + L + U + l)x{K+ K* + L + L* + 1) determinant with a block structure whose elements are the bi-orthogonal 
polynomials (pnrtpn ^nd their associated functions £,„, (see Il2.3b , (|2.4| |, (|2.11l l, ||2.12| |, (|2.13l l for their definitions) corresponding to 
the original weight. For example the bi-orthogonal polynomial 0„ is given by 



!=1 



(1.6) J](z-a,)n(l-«K')^«(^) 
= C-det ' 



>n-L'-L+c 



[a 



c=0,..,L-\ 



n-L'-L+cio^r)] r^i.-.K 

c=0,..,L-l 



[|S^-'^5„_L._L+,(|S,)] /=Y 



-0,..,L-1 

l,..,L 
-0,..,L-1 



=0,..,L-1 



[<pn-L'+c{z)]c=0,..,L'-'l 
[(pn-L'+c{ar)] 

[(p„-L'+c{a;)] 



t-a..,L*-l 

c=Q,..,L*-l 

I-I....L 
r=0,..,r-l 



[z- 



^n+c(z)]c=0,..,]iC* 
(>«+c(«r)] ' 



=1,..,K 
' c=0,..,K* 



[(^;)-5„«(^;)].^i 



a,..,K' 

^0,..,K* 



[Z (pn+K*+c(z)]c=l,..,K 

[(«;)- 



(>n+K-+c(ac)]j;;-'^ 
""^''t>n+X-+t:(«r)]7_\--«; 
[)S7^*5„+K.+c(iSr)].=Ci'' 

[(is;)-^^*^,,+K.+c(^;)W ^ 



The constant C is independent of z. The other components of the bi-orthogonal system are given by the above formula with the 
following substitutions: H/iiCz — II/LiCl ~ ctjz"^) 0*,(2^) is given by 1II.6I 1 with C 1-^ C*, 1-^ (\)\ in the first row, and second and 
third block rows and <J„ 1-^ <J* in the last two block rows; nf=i(z - f>i) n/=i(l - f>*-^~^) 3«(z) is given by (|1.6l l with c/)„(z) i-> <5„(z) in 
the first row; and nf=i(z - jS,) ny=i(l - /S]z~^) S* (z) is given by l|1.6l l with C 1-^ C*, (/)„ i-> in the first row, (\>n 1-^ <|)i^ in the second 
and third block rows and 4„ 1-^ 4* in the last two block rows. To find the leading coefficient K„ of the modified polynomials one 
can use the ratio formula 

.2 2[z«]0„(z) 2<1.;(0) 



[z«]H„(z) limz^ooH;,(z)' 



CJl 
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There is a recent parallel body of work jj), ||5l, l23l , | |24) which treats a BOPS 
defined with a measure d|U(x, y) on a two-dimensional domain (x, y) c and the 
rational modifications of this measure. This is more general in one sense than 
the problem we consider here and there are significant simplifications occurring 
in our system such as the existence of closed form expressions for the Christoffel- 
Darboux sums which are not present in theirs. As a consequence the resulting 
expressions have rather different forms. In the main these works only consider the 
normalisation constants of the modified system and are completely formal in their 
approach. In addition there is no attempt to develop a theory of the semi-classical 
weight case and its connection to isomonodromy preserving systems. 

For the class of regular semi-classical weights on the unit circle, taking for the 
purposes of illustration the simple example 



it was first shown in fl^ that this bi-orthogonal system is monodromy preserving 
in the spectral variable z with respect to deformations of the singularities zy. This 
fact is the unit circle analogue of the result for orthogonal polynomial systems 
on the line with regular semi-classical weights, first shown by Magnus in 1995 
p7J. Either of these systems constitute a realisation of a Garnier system and 
one can construct a dictionary to translate one into the other, and we refer the 
reader to [52 1 for details. In the regular semi-classical context the Christoffel- 
Geronimus-Uvarov transformations, under the specialisation of positioning the 
zeros or poles of the rational factors at the singular points of the weight, are the 
Schlesinger transformations of this Garnier system. The fundamental Schlesinger 
transformations shift the exponents pj pj±l, p* p*j^^r so that we are in effect 
considering a family of regular semi-classical weights with exponents displaced 
by integral values. 

The theory of Schlesinger transformations for Garnier systems was initiated by 
Schlesinger l43l and Garnier [18] where only the regular singular case was treated, 
and the theory for irregular singular points was completed in the contributions 
by the Kyoto school [28|,[29|,[30J,[33J,[32|. Schlesinger transformations are a spe- 
cial type of Backlimd transformation, which were systematically studied in the 
simplest non-trivial case of a Garnier system, i.e. the sixth Painleve system by 
Okamoto [40J. All the Backlund transformations can be composed of fundamen- 
tal transformations, generated by the elementary reflections and automorphisms 
of the extended affine Weyl groups oj/^for M = 3 and B^'^j for M > 4 and the 
permutation group Sm+i [36|,[27|,[49J,[47l. 

In the final section we present an alternative construction of the Schlesinger 
transformation theory to the studies mentioned above, starting from the BOPS 
with a regular semi-classical weight, which constitutes the most general classical 
solution to the Garnier system. Using all of the structures of such a system from the 
approximation theory point of view we can recover all of the known results, and 
some novel results to be reported in 1 52 j, in a very efficient and transparent manner. 
Furthermore the present study can also be seen as the completion of the tasks begun 
in liZJ, in that this work didn't cover the aspects of the integrable system relating 
to the Schlesinger transformations or the systems of bi-linear difference equations 
governing the T-functions. Our approach is to employ the simplest of classical 



(1.7) 
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arguments from approximation theory although we fully recognise that all of the 
results presented here could be found in an elegant manner using Riemann-Hilbert 
methods. A natural formulation of the bi-orthogonal system as a Riemann-Hilbert 
problem was given in [17| and it remains an interesting task to employ these 
techniques on the questions posed here. 

2. General Structures of Bi-orthogonality 

Let T denote the unit circle with C = e'^, 6 e (—1/ ti\. Consider an absolutely 
continuous weight w{z) and assume that its Fourier decomposition has meaning 
in a domain containing T, 

(2.1) w{z) = f wtz^ w, = f ^w{Q)Q-'. 

We define the bi-orthogonal polynomials {(/)„(z), (^„(z)}j^g with respect to the 
weight w{z) with support contained within the unit circle by the orthogonality 
relation 

(2.2) ^ _^j^(g,j;,,„(Q^„(g = 

Alternatively one can express this definition in terms of orthogonality with respect 
to the monomial basis 

(2.3) J^_|_.(C)^„(C)C«. 



(2.4) r ^w{QC"UQ = 






m <n 


1/k„ 


m = n 





m <n 


1/K„ 


m = n 



Notwithstanding the notation, cpn is not in general equal to the complex conjugate 
of (pn- We set 

(pn{z) = K„Z" + A„Z"-1 + llnZ"-^ + ... + (/)„(0), 

^„(z) = K„z" + A„z"-i + p„z"-2 + . . . + (^,,(0), 

where again A,„ p„, (/)«(0) are not in general equal to the corresponding complex 
conjugate. We have chosen the leading coefficients of both polynomials to be the 
same without any loss of generality as this preserves the relations arising from the 
orthogonal polynomial theory with, of course, the distinction made above. The 
initial values of the polynomials are c^o(z) = ^o(z) = kq and the normalisation 
implies WqK^ = 1. 

Denote the U{n) average by In[iv], 

n 

(2.5) I„lw] = {Y\ ^'^^^'^)u(n) " '^et[^''i-t];,it=o n-l- 

By a result of Baxter [4] the existence of the bi-orthogonal system is guaranteed if 
and only if I„[w] + for n e Z>o. We define a sequence of r-coefficients by 

(2.6) Tn = , r„ = , ro = rQ = 1, 
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which differs slightly from the standard definition of the reflection or Verblunsky 
coefficients a„ in that a„ = —fn+i ■ It is a well known result in the theory of Toeplitz 
determinants that 



(2.7) 



I„+i[w]I„-i[w] 



= 1 -rnrn- 



The coefficients are related by many coupled equations, two of the simplest being 



(2.8) 



kI = kI_^ + (p„{0)(pn{0), — - = r„r„-i. 



Introduce the reciprocal polynomial (p'^z) of the «th degree polynomial (p„{z) 

by 

(2.9) (/,;(z) := z>„(l/z). 
Fundamental to our study flTll is the matrix 

(2.10) V - IMZ) ^n{z)lw{z) 



(2.11) 



^"(^)--l</,;(z) -.£;(z)/Mz) ' 



where the associated functions are defined, for « > 0, 

C + 2 

1 



-wiQcpniQ), 



(2.12) ^n{z):=-z" f JLkLl^oiOcPniQ), 

Jt 



(2.13) 



T lni(, C - z 
rfC C + z 

27T/C C - Z 



A central object in the theory is the Caratheodory function, or generating fimction 
of the Toeplitz elements. 



(2.14) 



F(z) 



Jt 



C + z 



IniC, C - z 

The initial values of the associated functions are defined as <^o(z) = Kq[ivq + F(z)] 
and <J*(z) = Ko[wo - F(z)]. 

A consequence of the definitions (|2.3I (, ||2.4| | are the following determinantal 
representations for the polynomials. 



(2.15) 



(pniz) = — det 

in 



Wo 



Wn-1 
I 1 



W- 



Wn-j-l 

zi 



W-n\ 



W-i 

z" 



(2.16) 



( Wq ... W-„+i Z"\ 



'Aliz) = 7^ det 

-'n 



W, 



'<-! 



Wn 



W- 



y+i 



z' 



Wi I) 
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The associated functions have representations analogous to I l3.ll3.21 



(2.17) 



Uz) = ^det 



Wo 
Wn-l . . 

Igo(z) ■ • 

Wo 



W- 



W- 



Wn-j-l 



W-i 
gn{z)) 



(2.18) 



^;(,)._^det 



w 



n-j 



W„ 



W-n+1 gn{z) 



where 
(2.19) 



gj{z):=2z [ w{Q. 
Jt 



T 2niQ C - z 

A further consequence of the definitions is the difference system |19| 



(2.20) 

or its inverse 
(2.21) 

We also note that 
(2.22) 



1 / K„+iZ (pn+l{0) 



Y„^i :- K„Y„ - U^Q^^ ^^^^ 



Yn = — 



. Z Z 
-(p„+l{0) Kn+l ^ 



detY„(z) = - 



2z" 
w{z)' 



If one writes the coupled first order difference system ( |2.20t as one decoupled 
second order difference equation for the monic polynomial (pn '■- (pn/'Kn one finds 



(2.23) 



rn+i 

(pn+l (pn = Z 



(p„ (1 - r„r„)(p„-i 



These are the recurrence relations that serve to define one of the bi-orthogonal 
Laurent polynomials I34l , l35l , ll25l (the partner polynomial is foimd under the 
substitution r„ i-^ f„, (p„ (p„). 



Theorem 2.1 ( ||T9| ). The Casoratians of the solutions (p„, (p*^, E.„, are 
(2.24) 



^„+i(z)5„(z) - 6,+i(z)<|)„(z) = 2i^^^z«, 



(2.25) 
(2.26) 

valid for n>0. 



^„+i(0) 



.11+1 



(/,„(Z)C(Z) + UzWniz) = 2Z", 
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The bi-orthogonal system satisfies the spectral differential system |T7| (with the 
abbreviation d/dz := ') 



(2.27) -Yn-AnYn 

dz 



W(z) 



Q„(z) + V{z) - ^z0„(z) 



<^„+i(0) 



K„ 



(0) 



e„(z) 



ze;,(z) 



o;(z) - v{z) - ^@:{z) 



Kn Kn 

for n > under the sufficient conditions of the existence of the moments 



I 

Jt 



, [losro(Q]' - [logif(z)]' 

Q-z 



The particular parameterisation of A„ into the spectral coefficients W{z), V{z) and 
Q„(z), Q*(z), 0„(z),0*(z) will serve a purpose when we specialise to the regular 
semi-classical weights. We see that 



(2.28) 



. n w 
TrA„ = . 

z w 



The compatibility relation between the recurrence relations ||2.20| | and the spectral 
derivative (|2.27|l is 



(2.29) 



K'„ — A„+iK„ — K„A„, 



and this is equivalent to recurrences in n for the spectral coefficients, which can be 
foimd in [17J but are not required here. 

We will require the leading order terms in expansions of (pn{z), (()'„{z), <£„(z), ^* (z) 
both inside and outside the unit circle. 

Corollary 2.1. The bi-orthogonal polynomials (p„{z), (p*„{z) have the following expansions 
for \z\ < 1 



(2.30) (/)„(z) = c/)„(0) + (K„c/)„_i(0) + (/)„(0)A„-i)z 



Kn 



Kn-lK„-2 



{K„-l(pn-2{0) + (pn-l{0)A„-2) + 



Kn-1 



z^ + 0{z% 



(2.31) (p*„{z) = K„ + A„z + /i„z^ + O(z^), 

whilst the associated functions have the expansions 



(2.32) f .£„(z) = z" - —z' 

2 K„+l 



A„+lA„+2 lhi+2 



K»+lK„+2 K„+2 



,n+2 



An+lflji+S A„+3/i„+2 n„+3 A„+iA„+2A„ 



+3 



Kn+lK„+2 K„+2K„+3 K„+3 K„+iK„+2K„+3 



2«+3 _^ 0(Z"+*), 



+4\ 



(2.33) ^5;(z) = ^z"- - t!llim!!ll\,n.2 

2 K„+l \ K„+2 K„+iK„+2 / 

^n+3(0) _ (^n+l(0)p»+3 ^ An+3(0) / ^„+l(0)An+2 _ (^h+2(0) 
Kn+3 K„+lK„+3 K„+3 \ K„+lK„+2 K„+2 

T/je Zflr^e argument expansions \z\ > 1 ore 

(2.34) (pniz) = K„Z" + A„Z"-1 + finZ"-^ + 0(z"-3), 



2«+3 _^ o(z"+*). 



+4\ 



(2.35) (/);(z) = ^„(0)z" + ^(K„^„_i(0) + (^„(0)A„_i)z' 



n-1 



(k„-1(/)„-2(0) + (p„-l(0)A„-2) + 



K„-iK„-2 

whilst the associated functions have the expansions 



K„-l 



z"-2 + 0(z"-3), 



(2.36) 



(2.37) 



z + 



k2 (|)„+2(0) (/)„+l(0)A„ 



+1 



K^^l K„+2 



z-2 + 0(z-3), 



-«;(z) = 1 - 



Ak+1 _i 

z + 



An+2A„+1 



/in+2 \ 
K„+2 / 



z-2 + 0(z-3). 



\ K„+2K„+\ 

From the Riemann-Hilbert perspective [ 17| it is apparent that our bi-orthogonal 
system demands the existence of two domains, one containing z = Q and the other 
containing z = oo, where the above expansions are vahd in order to specify the 
system. 

3. Christoffel-Geronimus-Uvarov Transformations 

Starting with the determinantal representations of the polynomials and asso- 
ciated functions, given in l|2.15t , ll2.16b , (|2.17) l and (|2.18t , it is easy to recast these 
as multiple-integrals using the Heine formula. We observe that the following 
integrals occur with the weight modified by a rational factor 



(3.1) 



(3.2) 



(3.3) 



(3.4) 



/„[a;(Q(C -z)] = (-)"/„ MQ]- 



« > 0, 



I„[K;(Q(l-zri)]=I„[ii;(Q] 



</'«(z) 



Kfi 



n>0, 



InMQ 



I„[iv{Q) 



(1-zC-i) 
1 



r r /r\l ^ n I I ^ i 

lAMQ] — ^r:— i — / " > 0, |z| ?t i. 



] = (-)"/„ MQ] 



2z«-i 

K«-i'^;_i(z) 



n > 0,\z\ + 1. 



(C-z)' " ' ' " " 2z« 

These are the analogues of the Christoffel formula or the simplest cases of the 
Uvarov formulae for orthogonal polynomials |50|, |51|, and are specialisations of 
the Ismail and Ruedemann formulae for OPS on the imit circle |26| and those of 
Ruedemann [42 1 for BOPS on the unit circle. We seek to generalise the above 
expressions and derive formulae for the BOPS corresponding to an arbitrarily 
rationally modified weight l|1.5t in terms of the original BOPS. 

However Ismail and Ruedemann's study is incomplete for our purposes in a 
number of ways - only the modification of a real, positive weight on the unit circle 
by a quotient of a real, non-negative polynomial and a real, positive polynomial 
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is considered. This doesn't cover the case of complex weights, or equivalently 
non-hermitian Toeplitz determinants, and consequently the bi-orthogonal systems 
associated with this type of weight. In addition this type of transformation is a 
composite of fundamental transformations and it is not possible to dis-entangle 
these. In a subsequent work l42l the extension to the bi-orthogonal system was 
made, but still only treating the case of a composite transformation. In both of these 
studies formulae were given only for the modified pol5momials (p„, (p*„ whereas we 
require expressions for the complete system, that is to say for , as well. Another 
technical caveat applying to the results of f^S] is that the degree of the polynomial 
n cannot be less than the degree of the denominator polynomial L + L*. The studies 
by Godoy and Marcellan ll20l ,||2T | ,|22) are limited in similar ways, and have the 
additional complication that expressions for the modified orthogonal polynomial 
system are given in terms of determinants with elements containing reproduc- 
ing polynomials. Formulae for the polynomials and recurrence coefficients of the 
bi-orthogonal Laurent system imder a modification of the weight by multiplica- 
tion with a single factor, and division by a simple factor plus addition of a mass 
point were derived by Zhedanov Il54l . While we do not consider the problem of 
additional mass points at all here we will extend his results to arbitrary rational 
modification. 

We will carry out our task in seven steps, each one exhibiting our methods 
of deduction in their simplest setting, rather than give the proof for the general 
case. In the notation of the most general modification l|1.4|l, i ll. 51 we will derive 
formulae for the cases - K i= 0, n > L it n < L i= 0, K,L i= 0, K* i= 0, U i= 
and K*, L* + 0. The most general case has been given in 1II.6I 1. In each step we 
will denote the modified bi-orthogonal polynomials and associated functions by 
the same symbols <l'„(z), 0*(z), S„(z), H*(z) and it should be appreciated by the 
reader that the formulae only apply in their specific context. We also denote the 
modified leading and r-coefficients by K„,R„,R„; C, C* are arbitrary constants of 
proportionality independent of the spectral variable z; and lastly we denote the 
space of polynomials in z with degree at most n by n„(z). The same normalisation 
condition applies for the modified system, i.e. I|2.2b holds. We will assume hereafter 
that the following generic conditions apply - 

• that the zeros and poles are not pairwise coincident aj + a^, + a*. + 
a*^, |3* + jS*, ] + k,yn the proofs for convenience although confluent formulae 
can be derived from our results, 

• assume the existence of the original BOPS and of the modified BOPS, i.e. 
that r„f„ 7^ 0, 1 and R„R„ 7^ 0, 1 for all n e N. 



Proposition 3.1. For the polynomial modification of the weight 



(3.5) zv 



Yl(z- aj)zv{z), 



the corresponding bi-orthogonal polynomials are given hy 
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det 



( (p„{z) 

<t>n{0tl) 



(3.6) J](z - aj) (D„(z) = {-lfhn,„^K- 



det 



<t>n+K{z) 



,(pn{aK) ■ . ■ (p„+K{aK). 



(pn+K-l{ai)\ 



, n>0, 



.(pnioCK) ■ ■ ■ (pn+K-lioCK), 



det 



(3.7) Yl{z-aj)%{z) = {-l)%,„^K- 



det 



, n>0. 



and the associated functions 



(3.8) 



det 



( 5«(z) 



S„(Z) = (-l)Xn+K- 



det 



5n+x(z) 



,(/)„(ajc) . . . c/)„+x(ajc). 
'(]^„(ai) ... </)„+K-i(ai)' 



.(/'«(«k) • • • </>n+K-l(«K). 



, n > 0, 



(3.9) 



det 



s;(z) = (-i)^5„,„+K- 



det 



, n > 0. 
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Here b„^„+K = Kn/Kn+K, h„,n+K = ^n{0) / (pn+K{0) and the leading and r-coejficients are 
given by 



(3.10) 



and 



(3.11) 



and 



(3.12) 



det 



i=l 



det 



' K„ 


Kn+K ' 






4'„{aK) .. 




' 1/k„ .. 


1/k„+k ' 


(p'„iai) .. 




.(p'„{aK) ■ ■ 





det 



R„ = 



1 



(pn+K{ai) 



(p„{aK) ■■■ ipn+daK). 



1=1 «<■ 

det 



det 



K 

R« = (/)n+K(0)]^a,- 

!=1 



cpniai) 

(pn{aK) 



(pn+K-lioClY 
4>n+K-l{aK), 

</>;+x-i(«^)^ 



det 



Kn+K ^ 



Proof. Taking | |3.6|I first we note that nf=i(z - «/) ^n{z) e n„+x(z) so that 

K n+K 

(3.13) 



J](z - a^) cD„(z) = Vy<l>i (z), 

;=1 ;=0 



with 
(3.14) 



However by orthogonality | |2.4|I b„^j vanishes for j <n and in fact 

K n+K 

(3.15) J](z - Uj) 0„(z) = b„,j(Pj{z). 

j=l j=n 
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To determine the non-vanishing coefficients we observe that there are K vanishing 
conditions 



n+K 



(3.16) 



Y^b„,jcPj{ai) = 0, l = l,...,K, 



]=n 



which enable the {bn,j]"^^ ^ to be solved in terms of b„,„+K. This establishes 
Result (|3.8l l follows from (|3.6t upon using the definition (|2.11) l. 
To derive \3.7\ we start with 



n+K 



(3.17) 



j=l j=0 

For m = 0, . . . ,n — 1 we see that 



(3.18) 



(3.19) 







n+K ^ 



IniC 



-w 



C'^niQ, 



IniC. 



However the integral given above vanishes for j > n - m, and as 

(3.20) J _^^(QC(^^.(Q J _^^(QC"(^^.(Q, 

are non-zero we conclude that hn,j = for j - 0, . . . ,n — 1 and thus 



n+K 



(3.21) 



J](z-a,)cD;(z) = 2]^B„,,(/,;(z). 



]=n 



Again we have K relations 



n+K 



(3.22) 



Y^b„j(P'.{a,) = 0, l = l,...,K, 



j=n 



allowing us to solve for {b„_j)"^^ ^ in terms of bn^„+K- This yields jS.Vi . From this 
result we can derive l|3.9|l using the first definition l|2.12|l and | |2.9| |. This gives us 
the formula 



det 



(3.23) 



H;(z) = -{-lfbn,n+KZ"- 



cp:„{aK) 



but we rewrite the integral appearing here as 
(3.24) J 



dC C + z ; - - r dC P -z' - 



C-z 



+ Z' 
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n - zi 

The first term vanishes because (C + z) ^ _ ^ e ^jiQ and the second term evaluates 

to-z-«^;^/z). 

The relation for hn,n+K preceding ll3.10l l follows from examining the leading 
coefficient of Il3.6b and comparison with (|2.34l l. Similarly the relation for b„^n+K 
follows from the leading coefficient of l|3.7b and using | |2.35|I . The formula (|3.10l l 
is found from the leading coefficient of the z ^ oo expansion in ( |2.37t and the 
corresponding coefficient of the z ^ expansion in | |2.31| l. The relation ( 13. lit 
is derived from the ratio of the leading and trailing coefficients of | |3.6l l and the 
relation l|3.12b from the leading coefficient of the z — > oo expansion in l|2.37|l and the 
result for 6„^„+x- □ 

Proposition 3.2. For the reciprocal polynomial modification of the weight 



(3.25) 



10 



;z 



-w{z), 



n-=i(z-|S;) 

the corresponding bi-orthogonal polynomials are given by 



(3.26) 



det 



(-1)^ 

cI.„(z)=^^• 



det 



«+^-i£„_L(i3i) 



(| 3-"^^g„-,(|3,) ... rL"UPL)) 

1/k„ 



n>L, 



det 



(3.27) <l.;(z) = i-l^b' 



and the associated functions 



det 



det 



(3.28) J](z-^,)H„(z) = C 



' z^5„-l(z) 
|3-«+^£„_,(|Si) 



v/^l"^^'£"-l(^i.) 



n>L, 



|S-"^^5„_,(^i) ... |S-"-^5«-i(iSi)) 



, n>L, 



det 



)S-"-15„-i(^l)J 
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(3.29) 



det 



j](z-/5,)s:,(z) = (-i)X„z"- 



rz-"^%_^{z) ... z-"^;(2)^ 



, M > L. 



det 



Here b'„„ = Kn/K„ and the leading and r-coefficients are given by 



(3.30) 



and 



and 



(3.32) 



det 



K = (-)'- 



Kn-L 



^-«^^^„_l(jSi) 



i3r".£«(iSL), 



det 



1/K„ 1 



' 1 1 Kn-L 

^-"-^-15„_l(/^i) 

l,|3-"^^-l5„_L(iSL) ... rL~'Uh)} 



det 



(3.31) R„ = {-f<t>n{0)- 



^-«+^^„_L(iSi) ... ,S-«^i6,-i(|Si)^ 



det 



K„-L 



iS-""'««-L(iSi) 



t/5-«^^5„-L(/5L) 



det 



R„ = 



' (/)«-l(0) 



MO) ] 



K„ |'|3-"+^6,_l(jSi) ... f}-"^'U-i{lii)] 
det : : : 



Proof. We begin with the derivation of (|3.271 t. Writing 

II 

(3.33) 0,.iz) = Y^b:,j^jiz), 

j=0 
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we have 
(3.34) 



2niC 



1=1 



We observe that ni=i(C-j8;)</>/(C) e H^+lCO and therefore b;^. = OforO < / < n-l-L. 
Now consider the integral 

L 



(3.35) 







2niC 



-w 



for m - 1, . . . , L and subject to « > L. This can be rewritten 



(3.36) 
(3.37) 



" r 



j=n-L 



2niQ C - jS„ 



j=n-L 



which yields L conditions. These enable us to solve for V .in terms of 

n,n—L' ' n,n-l 

b'„ „ and we have the result Il3.27b . From this formula we proceed to the derivation 
of l|3.291 . Using the definition l|2.12|l on this last result we require an evaluation of 
the integral 



(3.38) 



r dc c + z 1 



We employ the partial fraction decomposition (assuming z i= jS;) 
C+z 1 C+z 1 1 



(3.39) 



^ 1 

T — 



1^-^"'(iS,„-z) nti^„(/Jm-iS;)' 



and evaluate the above integral as 



jSm + Z 



^iSr^"'C;(i3m)- 



When this result is inserted into the first row of the determinant formula | |3.27|I we 
see that the last term of this is a column independent linear combination of the last 
L rows of the determinant. Thus we can eliminate this term and arrive at l |3.29l l. 

To establish ( |3.26t we take an approach differing from the preceding arguments. 
We postulate the determinant 



(3.41) 



det 



'(pn-L+j{z) 

... |g-""'-^5„-L.;(iSi) .. 



-n+L- 



'5«-L^;(^l) 



19 



and examine its orthogonality properties with respect to the weight w 



<m <n and <k := L - j <L (we know n > k), and consider the integral 



. Let 



(3.42) 



-IV 



cv«-fc(c)e 



dC 



k-m 



n,=i(c-/?/) 



In the case m > k we have the partial fraction decomposition 



(3.43) 



£o 



^ 1 



N>0. 



c^n,=i(c-^o ' ^t-^nUp(^p-iS/) 

This can be inserted into the preceding integral and, by orthogonality and l|2.11|l , 
we see the first N + 1 terms vanish leaving 

-m+k-l ^ 



(3.44) 



p=l'nUp(|Sp-|S;)2 
In the remainder of cases, i.e. when < m <k, we use the decomposition 



(3.45) 



1 



Utilising this form in the integral we are left with 



N>0. 



(3.46) 



I 



k-m-l 



p=i'nUp(iSp-is,)2 



-<5„-t(jSp), n>m>k. 



In either case we conclude that integration of l l3.4Hl with the modified weight yields 
the determinant 



(3.47) 



det 



iS-«^'^5„_,(i3i) 



T^'^'^^n-kih) 



where the fp are independent of the column index. Therefore the sum in the first 
row is a linear combination of the lower L rows and the determinant vanishes. Thus 
0,1(2) is proportional to l|3.41t by ||2.3|I . All that remains is to settle the normalisation. 
To show that l|3.28t is true we utilise the foregoing result. Here we are lead to a 
determinant which has the following integral as entries of its first row. 



(3.48) 



^/c^„-k{Q. 



Then we can apply | |3.45| | to C'7 nf=i(C - §1) and so obtain the evaluation 



(3.49) 



L 



-/5p-«^*^5„-;c(/^p). 



By an identical argument to that employed with | |3.26|| we conclude that only the 
first term survives in the determinant and | |3.28|| follows. □ 
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We need to extend the preceding proposition to cover the case n < L, which is 
the subject of the following proposition. 

Proposition 3.3. For the reciprocal polynomial modification of the weight 

w{z) 

(3.50) w , -,2 = w 



ft 



n-=i(z-^;) 

with n < Lthe corresponding bi-orthogonal polynomials are given by 













z"<|)o(z) 




^h(z) 




(3.51) cD„(z) = C ■ det 






1 


«o(|Si) 


... /?- 








oL-n-1 

Pl 




1 
1 




•■• Pl 


'in(pL)) 






( ^ 







Z«(/)*(Z) 








(3.52) cd;,(z) = C • det 






^1 




... 
















■■• 






and the associated functions 




















' z^ 


-1 


... z« 


z"«o(z) 






L 

(3.53) [](z-/^,)S„(z) = 


C-det 




n-1 


... 1 


«o(iSi) 


... j3 










n-1 


... 1 

.. z«+i 


«o(|Sl) 


... fi 


«;(z) ) 


L 

(3.54) J](z-^,)S;(z) = 


C ■ det 




n 


.. |Si 




... ^ 








n 


.. |3l 


mo 


... ^ 





Proof. Orthogonality of the proposed formulae l|3.51t and | |3.52|| with respect to the 
modified weight according to the criteria l|2.3b and \2A\ is easily proved using the 
methods of the preceding proposition and the interpolation identity 

L gr 

(3.55) YL^^i—^ — =0, r = Q,...,L-l. 



p=l'n',p(^p-^;) 

The remaining formulae | |3.53| | and | |3.54| | then follow from these two again using 
the methods of the preceding proposition. □ 

Remark 3.1. There is a simple observation that renders the above formulae (13.511 - 
13.541 as transparently obvious - the inverse recurrence relation l|2.21l l allows for 
a backwards recurrence and therefore gives meaning to the bi-orthogonal system 
with negative indices n < 0. Recall that 



(3.56) 



Yn = 



( Kq{Wq + F) 

'^0 w 

Kq{wo - F) 



so that using l|2.21|l we find c^-i(z) = (p'_^{z) = and \k-iE,-\ = 1/z and \k-iE,*_^ = 1, 
assuming K-i i= 0. Proceeding in this way we next find that (p-2{z) ~ (p-i^^) ~ 
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and in fact all the polynomials with negative argument can be set to zero. The 
solutions of the backward recurrences for the associated functions can be shown 
to be pol5momials in 1 /z, namely with N e N 



(3.57) 
(3.58) 



1 ^ , , (p-N^m 1 
-k-n^-n(z) = + 



1 

rl^-N<^lA,(z) = 1 + • ■ • + 



^-N+l(0) 1 



yN-1' 



subject to the conditions k_jv, (p-^iO), (p-N{0) can be defined in some way and the 
K_jv are non-zero. We note that the above formulae are precisely the terminating 
forms of the general expansions Il2.32l2.33l2.36l2.37t . In defining the bi-orthogonal 
system with negative indices in this way we have to break the relationship between 
the polynomials and their associated functions as given in relations 112.1112. 1212. 13b . 
Then using elementary column operations on (|3. 51113.54) 1 we can recast these as 









„)(Z) 


zXz) ... 






(3.59) 0„(z) = 


C-det 




«)(iSi) .. 


«o(^i) ... 














^oiM ... 










( ^'^-(L- 




• z>o(2) ... 


(Pn(z) ' 




(3.60) <l);(z) = C* • det 


^l-"«*-(L 


















• mo •■■ 






and the associated fimctions 










I 

(3.61) J](z-j8y)H„ 


(z) 














( z^£.- 


(L-n)(z) 


. . . z"+i 


^_1(Z) z"5o(2) 


... Uz) ] 


= C • det 




-(L-«)(^l) 


... fil^ 


-i(iSi) «o(^i) 


■■■ /3-"6,(iSi) 






(l-«)(^l) 


... Pl^ 


-l(iSL) 5o(iSL) 


... f^rupii 



(3.62) J](z-|S,)H;(z) 



= C* • det 



^-(L-k)^ ' 



i^r"«i(L-„)(/^^) 



z"+ir_i(z) z"«o(z) 
^iSlidSi) 5*(iSi) 



r„(z) ^ 



Propositions 13.11 and 13 .21 can be combined into one result covering the case of 
rational modification of the weight. 
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Corollary 3.1. For a general rational modification, subject ton> L, 

njii(z-«y) 



(3.63) 



IV 



n -=1(2-13;) 



w{z), 



we have the corresponding bi-orthogonal system 
(3.64) 



ii[(z-a;)0„(z) = C-det 



(3.65) 



[](z-a;)0;(z) = C-det 



(3.67) 



a[(pn-L{ai) 


(pn{z) 

(pn{a\) 


(pn+K{z) ' 
(pn+K{ai) 


a^(pn-L{aK) 


(pn{aK) 


<pn+K{aK) 

iS^■£«+^c(iSl) 








«i</';_l(«i) ■■ 


(p;,{z) 


■ ^:.k(^) ] 




■ (pliaK) ■ ■ 






■ K'^l) ■■ 




' z^5«-l(z) 

«l(/)n-L(«l) 


(pn\(X\) 


5n+K(z) 
<pn+K{ai) 




(p„{aK) 


(pn+KioCK) 




■ ■ rLUh) 






■ (pniai) .. 






• (pliaK) •■ 

■ ^r«;(^i) ■• 











[](z-/?y)H;(z) = C*-det 



Proof. This follows by using the methods of Propositions 13.11 and 13.21 
We now present a "conjugated" version of Proposition l3.1l 



Proposition 3.4. For the polynomial modification of the weight 



(3.68) 



lU 



K' 



;2 



(3.69) 



J][(l-a*z-i)<l)„(z) = C-det 



(3.70) [|(1 - a*z-i) 0'„{z) = C ■ det 
and the associated functions 

(3.71) S„(z) = C • det 



K' 

]^(1 - a*z~^)w{z), 
/=i 




are ^z'uen by 






(p„{a\) .. 


z'^^'cp 


n+K-(z) 
/>n+K'(«P 








(m •• 

(/>;(«;) .. 













< «„(z) 



(3.72) 



S;(z) = C ■ det 



v'/'n(«K.) 



z-^^*^^„+K.(z) 1 

(a*)"*^*c/)„+K-(ai) 
(a;,,)-*^'c/)„+K-(a;^.)J 

(al)-^c/);,^.(a;) 



Proof. This result follows by employing the methods of Proposition 13. II 

There is also a "conjugated" analogue of Proposition 13.21 and is given in 
following statement. 

Proposition 3.5. For the reciprocal polynomial modification of the weight 

1 



(3.73) 



IV 



n^=i(i-is;z-i) 

the corresponding bi-orthogonal polynomials are given by 



(3.74) 0„(z) = C-det 



w{z), 



(3.75) f ;(z) = C ■ det 



det 


' i 


H-L-(Z) 

"5«-L.(|Sp . 


4 

■ rp- 


"5„(^P 






"««-L.(|SI.) . 






det 


( 


""Cl-^P ■ 


• • (isp 










•■ (iSI.) 





, n>V, 
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and the associated functions 



(3.76) fj(l-^;z-i)S„(z) = C-det 



(3.77) 

L* 



J](l-)S*z-i)H;(z) = C-det 



^ «„-L.(z) 

(^;)-"5„-L.(^;) 

Ui3-.)-"5„-L.(iSL.) 



5«(z) ■ 

(is;)-"5„(is;) 

(isi.)-"5„(isi.), 

r„(z) 
(/3;)-"5;(^;) 



n > V, 



n > V. 



Proof. This result follows by employing the methods of Proposition l3.2l □ 

Analogous to Corollary 13. 11 Propositions 13. 41 and l3.5l can be combined into one 
result covering the case of a "conjugated" rational modification of the weight. 

Corollary 3.2. For a general rational modification, subject ton> L* , 

nf:i(i-a;z-i) 



w 



K' 



(3.78) 

we have the new bi-orthogonal system 

K' 

(3.79) J](l-a;.z-i)cl.„(z) 



(pn-L-{z) 
(pn-L'{a{) 



n^:i(i-^;z-i) 



w{z), 



C-det 



(^;)-«6,-L.(^l) 



(p„{a{) 



Z '^'(pn+K-iz) 

(ap"*^*(/)„+K.(a*) 



U|S*.)-«6,-L.(^I.) ... ifil.rUPl.) ... (^'J-'-'^'WdSIJj 



K' 



(3.80) [](l-a;.z-i)<l.;,(z) 



= C ■ det 



(/);(z) 
(/);,(«;) 



(|SI.)-"Cl-(^I-) ■•■ (iSI.)-".£;,(i3I.) ... 
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(3.81) fj(l-/^;z-i)H„(z) 

(pn-L-{a{) 



= C ■ det 



(pn-L-(a'i..) 

(/5;)-«6,-L.(/31) 



(a^.)"^>«+K-(«K-) 

(,s;)-«-^^'5„^K.(|s;) 



(3.82) fj(l-)3;.2-i)H;(z) 



= C* . det 



(/3i)-".£:,_,.(/3i) 



(4.)"*'>;+K.(4.) 



Proof. These results may be derived by utilising the methods of Propositions 13.11 
andlMl □ 



Corollaries l3.1l and l3.2l can be combined into a single result for the most general 
rational modification of the weight (|1.51 and this is stated in (|1.6|l without proof. 



Remark 3.2. We now show how the results of Ismail and Ruedemarm 11261 and of 
Ruedemann 1 42 1 can be recovered from a specialisation of the result for K, K* > 0, 
which follows easily by combining | |3.6|I and l|3.69b . 



(3.83) [](z-a,)J](z-a;.)0„(z) 

;=1 ;=1 



C-det 



{aKf (j)„{aK) 
(a;)^>„(a;) 



(pn+K'{z) 


■ (pn+K'+K{z) ' 


(pn+K'{ai) ■■ 


. (pn+K'+K{ai) 


4>n+K'{aK) ■ ■ 


■ (pn+K'+KioiK) 


<t>n+K'{a\) .. 


■ (l)n+K'+K{a{) 







, n>0. 
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The other result of \26l is a likewise a special case of L, L* > 



(3.84) cl)„(z) 



= C-det 



(is;)-""^6,-L-L.(^;) 



U^2.)-'""««-L-L.(iSI.) 



(/>n-L-(z) 

(^l)-«5„-L.(^l) .. 

(^;)-"£„_,.(^;) .. 

(^I.)-"5„-L.(^I.) .. 



(p„{z) 
(^i)-"5„(iSi) 



, n > L+V 



In the context of orthogonal polynomials Theorem 1 of ||26l states a formula which 
is extended to the bi-orthogonal polynomial setting by Theorem 2 of [42 1. This 
is the K - K* = m case in (|3.83t and we now proceed to show that our result is 
equivalent to theirs. Our result can be expressed, keeping only the first row, as 



det(z"'(/)„(z) zcp (Z) <P n+2m 



and we use elementary column operations and the forward and backward recur- 
rence relations l|2.20b and l l2.2Hl . The initial step is to recast this as 

det(z'" ^(p*„+i, ■ • ■ , (pn+m' (pn+m, (pn+m+1, (pn+2m), 

and then conduct a sequence of m consecutive sweeps. The first sweep uses (p*„+,„ 
and those on and to the right of (pn+m+i in a sequence of successive steps. At 
each of these steps a recurrence relation is used to transform the first element by 
incrementing its index and then another is used to decrement the index of the 
second and thereby introducing a factor of z. (we use • to denote an intervening 
element not involved in the column operations) 

det(z'" ^<p*„+i, (pn+m' '' <pn+m+l, <pn+2m), 
i 

det(z'""^(/)* . . . , (/)*+,„+!, •, (pn+m+l, (pn+2m), 
i 

det(z'"-l(/);^^, . . . , ■, Z(pn+m, (p„+2m), 

i 

det(z'" "^<^*+i, . ■ . , (pn+m+2' ■' ^'Pn+m, (pn+m+2, ■■■ , <pn+2m), 

i 

det(z™ . . ■ , (pn+m+2' ■' ^'Pn+m, ^'pn+m+l, ■ ■, <pn+2ni), 

det(z"'"^(^*^j, . . . , (/)*+2m-l' ■' ^'Pn+m, • ■ ■ / Z(p„+2m-l)- 



27 



Then we successively use a recurrence involving the first element and the second 
in reverse direction in order to decrement its index, restoring the index ion + m. 

det(z'" ^<p*,i+i, ■■■ , (pn+2m-l' ^n+nif ZCpn+m, Z(pn+lm-\)i 



det(z"' • • • ' ^n+m' (pn+m, Z(pn+m, Z(pn+2m-l), 

The second sweep uses zc/)* and those of z(p„+m+i and the ones to the right in a 
similar sequence of steps 

det(z'«-V;+i z(p Z(pn+m+l, • • • / Z(()„+2m-l), 

I 

det(z'"" (/)* . . . , z(p*„+„, ■, ■, z(p„+m, 

i 

det(z'""^(/)*^j, . . . ,Z(^*+^, •, ■,Z(p„+m,Z^(p„+m, ■ ■ . , Z(p„+2m-l), 



det(z'" ^(p„+ir ■ ■ ■ I Z(pn+2m-2' "' '^'Pn+mi Z,^(pn+m/ ^^n+m+li ■ • ■ i Z^</>n+2m-2)/ 
detCz^'^e^* . . . , Z^;+„, •, •, Z<^„+m, Z^(^„+m, . . . , Z^(/)„+2m-2), 

We repeat this procedure m times and the final sweep yields 

det(z'" ^(pn+m' ■ ■ ■ I Z'Pn+m' 'Pn+m' 'pn+miZ<pn+miZ'^(pn+mi ■ ■ ■ /Z'"(pn+m)/ 

which is the desired expression. 

We now specialise the above results and consider the K = 1,L = and K = 0,L = 
1 cases in further detail as these particular transformations form the generators 
for the group of Christoffel-Geronimus-Uvarov transformations. To this purpose 
we define the matrix generators R^{z;a) and K~(z;j3) for the K = 1,L = and 
K = 0, L = 1 cases respectively 



(3.85) y:(z) = Y„ 

(3.86) Y-{z) = Y„ 



1 

;z 







1'" 



= Kiz;a)Y„{zl 
= R;(z;i3)y„(z). 



In our specialisation to the regular semi-classical weights we will draw heavily 
on the results for the K = 1,L = and K = 0,L = 1 cases as they will form the 
generators of the Schlesinger transformations in that context. 

Corollary 3.3. For K = 1,L = with the numerator factor ofz - a denote the coefficients 
of the transformed BOPS by K'l,,r^,r^ and the bi-orthogonal polynomials and associated 
functions by <!?„, O*, H„, H* . Assuming (pn{a) + Ofor n e Z>o these are given in terms of 
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the base system by 

(3.87) 

(3.88) 

(3.89) 



+ \2 



r — 



^ (p„{a) 
'^/j+iKfj— - -, 
9n+i(a) 

niO)(pn+l{a) - (pn+l{0)(pn{a) 



(p„{a)' 



and the corresponding bi-orthogonal polynomials and associated functions are 



(3.90) 0„(z) 

(3.91) o;(z) 

(3.92) S„(z) 

(3.93) S;(z) 



(p„{z) ■ 



(|)„+i(0) c/);,(a)(/)H(z) -(/)„(«)(/); (z) 



K„+i(p„{a) 



z — a 



z - a 



K„ 

1'+ 

K„ 



a </'„(«)</'« (z) -</'«(«)(^„(z) 

——-(/)« (z) + 

(/>n(a) 

(z-a)<5„(z) 



(pn{a) 

K„+l(pnia) 



{z-a)^„{z)- 



(p„{a) (p„{a) 
Alternatively the matrix generator is given by 



[(p'MKniz) + (/)„(a)r„(z)) 



(3.94) 



K{z;a) 



K„ z - a 



( _ Kn(p„+l{a) (pn+\iO) \ 

K„+i(|)„(a) Kn+l 

^Z 

(pn{a) 



Proof. This follows easily from Proposition 13.11 and the expansion formulae l|2.30b - 
(IZ37t . □ 

Remark 3.3. This case has been given by Zhedanov fSl] in the context of BOLPS, 
which he termed as a Christoffel transformation, so we can directly compare some 
of the above results with his. Accordingly we verify that Eq. (2.6) and (2.7) of l,54J 
follow from lll87t and lIlM , while Eq. (2.1) follows from JS^. 

Corollary 3.4. ForK = 0,L = 1 with the denominator factor of z- f} denote the coefficients 
of the transformed BOPS by K~,r~,r~. Assuming £,*„((}) + Qfor n e Z>o these are given 
by 



(3.95) 
(3.96) 

(3.97) 



(K„) — KnKn-l ^* (^g) ' 



r = 



r = 
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and the corresponding bi-orthogonal polynomials and associated functions are 
(3.98) 



0«(z) = 
(3.99) 

3>;(z) = 
(3.100) 

H„(z) = 
(3.101) 

2;(z) 



[e„{ii)u^) + uiiWn{^)]. 



K„ 
K„ 



^ (z - j3) ( ^c/)„(z) - c/,;(z) I + (5;(iS)(/)„(z) + urnii^)) 



Kn p .^;,(is)«„(z) - 6,(^)«;(z) 



Z-fi 



5„(Z) + ^5„(Z) + ^-^ 



Alternatively the matrix generator R„ is given by 



(3.102) 



^■£;(|S) 



Proof. This also follows from Proposition 13.21 and the expansion formulae (|2.30l l- 
(I237t . □ 

Remark 3.4. This case is specialisation of the Geronimus transformation given by 
Zhedanov |54!| in the context of BOLPS. These transformations include a contribu- 
tion from a mass point which we do not consider here, and we effect this by the 
choice of the variable (in his notation) /,( = |S, (p„ = l;„{fi)/Kn. Consequently we find 
V„ = /3Kng„-i()g)/ (|)„(0) g;_i(j 3) and therefore we verify that Eq. ( 2.15) and (2.16) of 
11531 follow from ( |T95l and i33B , while Eq. (2.11) follows from | |198l l. 



Compatibility of the Christoffel-Geronimus-Uvarov transformations (|3. 8513. 86b 
with the other fundamental structures of the bi-orthogonal system imply a number 
of simple relations. For example compatibility with the n n + 1 recurrence (|2.20t 
implies the following result. 

Proposition 3.6. The K„ and R„ matrices satisfy the compatibility relation 
(3.103) R„^i(z)X„(z) = ^C,t(z)R^(z), 

for all z, n, a,li. 

Proof. Considering the X = 1,L = case first we use the explicit forms for the 
matrices K„,K'^, given by l|2.20l l along with | |3.87I3.88|I , and the matrices KJ,RjJ"^^ 
as given in (|3.94t . A direct calculation of the difference of the left-hand and 
right-hand sides then yields the null matrix upon repeated use of the recurrence 
relations for (pn+iic^), (pn+i{a), (pn+ii^) implied by | |2.20| |. For the X = 0, L = 1 case we 
used llZ20l l and Il3.95l3.96b for the matrix K' and the formula (|3.102b for R;,,R~^y 
The appropriate difference was shown to be identically zero using the forward 
recurrence for 4* ( jg) as implied by ( 12. 20b and the backward recurrence for 4* _j (|S) 
as implied by l|2.2Hl . □ 
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The other structure that we have to verify compatibihty with the Christoffel- 
Geronimus-Uvarov transformations is the spectral derivative, namely (|2.27) l. In 
order to verify this we have to find expressions for the transformation of the spectral 
matrices ^^^(z) under the action of the generators X = 1,L = and K = 0,L = 1, 
which are given in the following lemma. 



Lemma 3.1. With the assumptions of Proposition 13.31 the transformed spectral matrix 
A'^{z) under the Christojfel-Geronimus-Uvarov transformation K = 1,L = is given by 

(3.104) (z-a)A:(z) = — 1— f-'^-^W"^'^^ MOma) 
" K„(p„+i{a) \ aK„+i(p„(a) -axn+icpnia) 

%,+iia) - K„+i(p„{a)z -(p„+i{0)(p„{a)\^ 



K„K„+i(p„{a)(p„+i{a) \ -Kn+i(p*„{a)z aKn+i(p„{a) 

-aK„+i(p„{a) -(pn+\{0)(pn{a) \ 
-Kn+i(p'n{a)z K„+i(p„{a)z - K„<p„+i{a)j ' 



or 



(3.105) A:(z) = + K{z)A,{z)[R:iz)]-\ 

and the analogous result for A~{z) under the Christoffel-Geronimus-Uvarov transformation 
K = 0,L = 1 assuming the conditions of Proposition \3.4\ is given by 

(3.106) iz-p)A-M = m^,^ I 



jSk„-ik„.£;_i(jS).£*(jS) \(|)„(0)5;(^5)z -K„5;,(/^)z + ^k„_i5;_i(/5) 



K„5;(,s)z - /3k„-i5;_i(/^) ^^nUI^) 

^' <^«(0)5;(|S)Z -liKn^M' 



or 



+ R,;(z)A„(z)[R,;(z)]-i. 

z=/i 



Proof. To establish the first result | |3.104| | we differentiate the expressions ( I3.87|3.88t 
with respect to z and employ the matrix equation i2.27} for the derivatives of (p*„ . 
Having rendered the two derivatives in terms of (p,, and (p*^ we invert l|3.87b and 
(|3.88l l so that these two polynomials are given in terms of 0„ and O* . The second 
result is found by an identical argument starting with (|3.95l l and ll3.96l l instead. □ 

Now compatibility of the Christoffel-Geronimus-Uvarov transformations with 
the spectral derivative must imply the following result. 

Proposition 3.7. Compatibility of the Christoffel-Geronimus-Uvarov transformations 
K = 1,L = and K = 0,L = 1 with the spectral derivative imply the condition 

(3.108) ^Rtiz) + Rt{z)A„{z) = At{z)Rt{z), 

for all n, z, a,fi. 
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Proof. We shall verify (|3.108|l by treating the ± cases separately. For the + case 
we employ (|3.94t , which we also differentiate with respect to z, and the formula 
for A'I'j ||3.104| |. Then it is a matter of an exercise to check that the difference 
between the left-hand and right-hand sides are identically zero without making 
any assumptions on A„. The — case is treated in a similar manner starting with 
(|3l02l l and (|3l06l l. □ 

In concluding this section we shall display the matrix generators for modifi- 
cations with X* = 1 or L* = 1 for the sake of completeness even though these can 
be found from composite transformations with K - 1,L - 1 and [} = 0, a - 
respectively These results also exhibit an elegant symmetry with respect to the 
other generators. 

Corollary 3.5. For X* = 1, L* = with the numerator factor of 1 - a*z~^ denote the 
coefficients of the transformed BOPS by k^, r^, and the bi-orthogonal polynomials and 
associated functions by 0„, O* , H„, S* . Assuming (p*„{a*) + Ofor n e Z>o these are given 
in terms of the base system by 



(3.109) 
(3.110) 
(3.111) 



(k„) = Kn+\Kn 



r — 



(pnia'Y 



and the corresponding matrix generator is given by 

( 



(3.112) 



<(z;a) = 



K„z-a* (p„+i{0)a' ^ K„(/);^i(a*) __ ^ 

K„+i (/);(«*) 

Proof This follows from Proposition |33| and the expansion formulae l|2.30|l - (|2.37l l. 

□ 

Corollary 3.6. For K* = Q,L* = 1 with the denominator factor ofl- fi*z~^ denote the 
coefficients of the transformed BOPS by k~, r~, f~. Assuming ^n{fi*) + Ofor n e Z>o these 
are given by 



(3.113) 
(3.114) 
(3.115) 



(k„ ) — KnKn 



K„<£„-l(/S*) 



5;-i(r) 



|3*5«-i(iS*)' 

and the corresponding matrix generator is given by 

K„-l<£„-l(jS*) 

^'»r„(r) 



(3.116) 



K„-l£,„-l{p') 

y K„_ijS*5„_i(|3*) K„_i|S*5„-i()S*) 
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Proof. This follows from Proposition |33] and the expansion formulae l|2.30b - (|2.37|l . 

□ 

4. The Semi-classical Class of Weights 

From the viewpoint of classical orthogonal polynomials the weight (11.31 1 is 
a particular example of the regular semi-classical class | |4.3K characterised by a 
special structure of their logarithmic derivatives 

1 d 2V(z) Pi 

Here y(z), W(z) are polynomials with deg V{z) < M+l,deg W{z) = M+1. The zeros 
of W{z) define M movable singularities |zy}^^ located at arbitrary positions in the 

finite complex plane. The residues at the finite singularities, the parameters {p/l^g, 
are also arbitrary complex constants only restricted to ensure that the trigonometric 
moments have meaning. In addition the matrix spectral differential equation in z 
(|2.27b has a distinct set of fixed singularities, due to the bi-orthogonal system and 
its support on the unit circle, at Zq = and Zm_+i - °° for generic values of the 
parameters pj. It is possible that one of the zeros of W{z) coincides with the origin 
and we cover this situation by denoting this zero by zq and its exponent by pg. 

To avoid technical complications we require the following generic conditions 
for the regular semi-classical class - 

(1) deg (W) > 2, 

(2) deg (V) < deg (W), 

(3) the M + 1 zeros of W(z), {zo,zi, . . .,zm] are distinct, 

(4) the residues pj = 2V{zj)/W'{zj) ^ Z>o. 

From the expansion of the denominator in terms of elementary symmetric functions 

M M+1 

(4.2) W(z) = - ^/(*)) = £(-l)'e;[z,]z*^^i-'. 

j=0 1=0 

and (|4.1|l one sees our system is characterised by the singularity data {zy, pyj^p 
of the weight, and our task is to deduce the consequences for the bi-orthogonal 
system. If | |4.1| | is viewed as an ordinary differential equation for w{z) given the 
data {Zj, p/I^Q then one particular solution is the generalised Jacobi weight 

M 

(4.3) Mz) = J](z-z//, pjeC. 

j=Q 

We will find that a certain symmetry will be present if we consider the following 
"conjugate" expression to the above weight 

M 

(4.4) w{z) = Yl{l-ZjZ-Y', P*eC. 

/=i 

For the regular semi-classical class the spectral coefficients 0n(z), 0*,(z), Q„(z) 
and O* (z) in (|2.27|l are polynomials in z of fixed degree deg Q„ (z) = deg O* (z) = M, 
deg 0„(z) = deg 0*,(z) = M - 1, independent of n. This is true provided that 
W(0) = 0, i.e. one of the singularities is located at the origin zq = 0, and is taken 
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only in order to simplify the ensuing expressions. Once the singularity data is given 
these polynomials themselves are most usefully characterised in two distinct ways 
- either in terms of the coefficients of the monomials in z or their evaluations at the 
singular points {&„{zj), &n{zj), Q„(z^), D*„{zj)]^^. 

As a consequence of the polynomial form of the 0„, 0*,Q„, O* the spectral 
matrix A„ appearing in l|2.27|l has a partial fraction decomposition 



M 



(4.5) 



j=0 



Z - Z; 



in terms of the residue matrices A„j. The «-th residue matrix A„j at the finite 
singularity Zj is given by 



An,j — 



Pj 



2V{zj) 



Kn+1 



-Q„{Zj) - V{Zj) + Zye„(Z;) 



(^«+l(0) 



Z;0;,(Z;) 



-e„(zy) 
q;(z;) - v(zj) - ^o;(zy) 



(Af.-) ^ _fL( (PniZi)^n{Zj) -(p„{Zj)^„{Zj)\ 
^^■"^ 2z^\-cP'„{Zj)^'„{Zj) (P„{Zj)i'„{Zj)j' 

for j = 1, . . . ,M, while for / = the expression is 

(4.7) A„,o = (n-po)(j 7' 
and for the singular point zm+i = ^ 

(4.8) A„,oo 







-n 



Let us define the diagonal matrices of formal monodromy Ty for = 0, 1, . . . , M, oo 
by the diagonalisation 

(4.9) An,j = Gn,jTj(Gn,i)~\ 

and we compute these as 

To ^ 

[o n-pa' 

j = l,...,M, 



(4.10) 



7 = 




-Pi) 
-n 

M 



] 



(4.11) 



Cn,j — 



II T^oP,} 

(pn{0) dg/Kn 

Kn ' 

(f)n{Zj) djE,n(Zj) ' 

<Pl{Zj) -dj^liZj)^ 

Kn 

j^^„(0) doo/jC,,^ 



7 = 
, i = l,...,M, 

j = oo 
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where do, di,. .., du, doa are fixed by the weight and its support. The reader should 
note that our conventions are sHghtly different from those of the Kyoto School in 
that detA„y vanish for j - 0, . . .,M and TrA„ j are non-zero in the generic case for 
i = 0, ...,M, oo. 

There are sum identities implied by the above relations (|4.6|l and 1 14.81 



Pk = {n-po)r, 

pk 2? " 

k=l k 

(p„(Zk)C„{Zk) 
pk 

k=l k k=Q 

'PniZkWniZk) [ V 1- 
Pk = \^ + l^Pk r, 

k=l k \ k=Q ) 



What is clear from our experience is that the evaluations of the polynomials 
&„, 0* , Q„, Q* at the singular points and their subsequent use as auxiliary variables 
greatly simplifies the theory. We shall therefore seek to cast our theory in terms 
of these quantities. Firstly we note that the bilinear products of a polynomial and 
an associated function evaluated at the finite singular point Zj are related to the 
spectral coefficient evaluations at the same point by the following equations, valid 
for all; = 1,...,M 



(4.16) 
(4.17) 



rf)„+i(0) 0„(z,) 

Uz,,4.(z,) = 2*-ii.;^. 



Kn ' 2V{Z jY 
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(4.18) 
(4.19) 
(4.20) 
(4.21) 
(4.22) 
(4.23) 
(4.24) 
(4.25) 



(h„+l(0) ^, □;(Z;) + V{Zi) 

(b*(z)r (z) - 2 2"+^ — 

'P„lZ;)4„,ilZ;) ^ Z^- 2y(z^.) ' 

4)n+l{0) ^.□;(Z/)- V(Z,) 

c/);,,(z,.).£;,(z,.) = -2ii^z"-i "^'^ 



(/)„(Zy).S;,(Zy) 



V(Z;) L 
Z 



2V{zj) 

0„izj) - V{zj) - ^^Zj&„{zj) 

Kf, 



1 



V{zj) L 
z" 



o;(z,)-y(z,)-^e;(z,) 



V{zj) L 

z' 

V(Z;) L 



CLniZj) + y(z^) - ^Zy0„(Zy) 



q;(z,) + y(z,) - ^e;(z,) 



In the semi-classical class the singular points are free parameters in the com- 
plex plane and another key structure in our integrable system are the deformation 
derivatives of the system | |2.10| | with respect to arbitrary trajectories of the sin- 
gularities Zj{t),j = 1,...,M parameterised by the deformation variable t. The 
deformation derivatives were derived in I.17J , and given by (where d/dt := ') 



(4.26) 



y - R Y 

at 



M 

L 



A 



Z - Z; 



with 



(4.27) 



Boo — 



Kn 







K„(p„{0) + Kn(pn{0) _Kn 



This latter term arises in our system because the leading term of Y„ as z — > oo is 
not normalised to the identity matrix as is the case in most other treatments. 

Compatibility between the set of deformation derivatives and the spectral 
derivatives l|2.27b leads to the system of integrable non-linear partial differential 
equations for the residue matrices known as the Schlesinger equations. 



(4.28) 
(4.29) 



An,j = [Boo, A„,j\ + ^^^—^ [An,k, An,j] , j = l,...,M, 
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A number of deformation derivatives v^ill be required subsequently, which v^ere 
computed in (TT). 



'Kn f 

h 

M 

(4.31) '"=Lr^., 2z" 



(4.30) 2^ = -Z 



PK 



— Zk 2z« 

k=l k 

PK 

M 



(4.32) ^r,=Yj Pk 



Zk [r„(\)„{zk) - (\)*„{zk)\£;,,{zk) 



— Zk 2z" 
;c=i 

In addition we will need the deformation derivatives of the polynomials and 
associated functions evaluated at a singular point P„y := (p*„(zj)/(()„{zj) and Q„j := 

UZiV^riZj) 



(4.33) Pn.i^Yu 



Pk 



Zk [(pn{^k)Pn,j - (p'„{ZkWn{Zk) 



Zk 2z" 

k=l k 



Zj-Zk [ct)„{Zk)P„,j - (p*„{Zk)][£.n{Zk)Pn,j + 5*,(Z/c)] 



and 



(4.34) Q„,y = -£ 



-, Zk 2z" 

k^j=0 I " k 



^ Zk QnMni^k) + 't>'n{Zk)Qn,i\ini^k) 

Zk 2z« 



Pk- 



k=l ^'^k 



pk: 



Zj - Zk [(pniZk) + <p;,{Zk)Qn,j][£.n{Zk) - 5;,(Zjt)Q„,;] 



Zk 2z« 

k*j=0 ' k 



which can easily be deduced from the above work. 

A consequence of the evolution of the residue matrices under ( 14.281 1 is the 
result that arbitrary deformations preserve the monodromy data, in this case the 
collection of monodromy matrices {Mj]^^, with each defined by 



(4.35) y„ {zj + 6e^"') = Y„ (zy + 6)My. 

Each monodromy matrix has an upper triangular form, one of the three structures 
expected for classical solutions of the Gamier systems l38l , 

(4.36) M, = C-V-^'C, = (J 7 = M, 

where the cy, and thus the My, are independent of Zk or t, and also of n. The 
connection matrices are given by 

(4.37) C; = |J '/J, / = 0,...,M. 

We note that the BOPS satisfies the necessary condition for being a classical so- 
lution of a Garnier system 1.52 J . This means that the formal monodromy exponents 



37 



given by 

(4.38) 00 = n - po, 

(4.39) Oj^-pj, 7 = 1,...,M, 

M 

(4.40) e^^n + J^ pj, 

satisfy the classical condition Yjo,i m,cx) = 

There is a final structure in our rntegrable system that needs to be put into 
place and this concerns the discrete deformations of the parameter set of the sin- 
gularity data, the set of exponents (p/l^g. These are known as the Schlesinger 
transformations and are the subject of the next section. 



5. The Schlesinger Transformations 

When one of the free parameters a, jS, a*, jS* in Propositions 13.31 13.41 13.51 or 13.61 
co-incides with a singularity of a regular semi-classical weight, say zy, then the 
weight l|4.3t undergoes a modification pj pj ±1 or p* i-^ p* ± 1 respectively. In 
this situation we are in reality treating a family of weights whose exponents differ 
by integers and in the context of isomonodromic systems these modifications are 
Schlesinger transformations. These transformations and the basic recurrence entail 
the following types of transformations in terms of monodromy exponents 



n 


n±l 


do 1-^ f 


%±1A 




3oo 


± 1 


Pi^ 


pj±l 


dj 1-^ f 


3y + 1 , 6 






+ 1 




p;±i 


00 ^ 


00 + 1, 


dj 




Fl 


n n±l , 


Pj ^pj + 1 


00 ^ 


00 + 1, 


dj 




tl 


n n±l , 


p*^ p* + l 


dj 1-^ f 


3y ± 1 , 6 




9oo 


+ 1 



If we define the matrix transformations for the mappings py i-^ py + 1, p* i-^ p*±l 
for; = 0,...,M, 

(5.1) y„(z; Pj ± 1) = Rr(z; P;)y„(z; pj), 

(5.2) y„(z; p* ± 1) = Rf (z; p,)Y„(z; p*), 

(cf . I l3.85l l,l l3l86t ) then these are given in the following proposition. For this and the 
following propositions it is only necessary to exhibit proofs for the cases as the 
R*J~ transformations are composed from the R'^ and RjJ'^ transformations. 

Proposition 5.1. The elementary Schlesinger transformations for the shift pj ^ pj ± 1, 
with Zj + 0, are given in terms of the base system with parameter exponents {pyl'^g by the 
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formulae 
(5.3) 

(5.4) 



1 < 

2 2 1 )Cj^ 



( _ KnCpn+ljZj) (/)„+! (0) ^ 
^ K„+i(/)„(Zy) K„+i 

-6 ; 



K„-i^;_i(z;) K„-i<5;_i(z^) 

'?'«(0)r„(Z;) 2 K„5;(Z;) 2 



U«-l5^-l(Z;) ^/ K„-i^;_i(Zy) Z; j 

assuming that (p„{zj) + and ^*„{zj) + Qfor all n e Z>o respectively. 

Proof. The results (|5.3t and ||5.4| | follow immediately from (|3.94l l and l l3.102t respec- 
tively. □ 

Remark 5.1. One can verify 



(5.5) 



detRr(z;P;) 



z - z 



■, detR/, (z;p;) = z-zy. 



using the forward recurrence for (p„+i{zj), i.e. | |2.20t , in the first equality and in 
the case of the second, the forward recurrence for ■5n(zy) | |2.20| | and the backward 
recurrence for E„-i{Zj) l|2.21b plus i2.8[ . 

Corollary 5.1. The transformed leading coefficients are given now by 



(5.6) {xlf = -K„+iK„ 



(pn+l(Zj) 



e„(zy) 



Q„(z^) + y(z;)' 



ZjKn+\(pn{Zj) K„+i 



(5.7) rl = 

(5.8) = 

and 

7 JZi) 0* ,(z,) 

(5.9) (k„) = -K„K„-1— ;— = -K„K„-iZ^ — 



Zy0«(Z;) 



^«(Z;) K„ Q. (2^.) _ y(2^.) _ ^e;(z,) 



(5.10) 
(5.11) 



r = 



-£;,(z;) 

zy<£„-i(Z;) (/)„(0) 



n;_l(Z;) + V(z,)' 
Zye„-i(Zy) 



(Z;) K„_i Q„_j (2 .) - y(2^.) - j^Z^0„_i (z^) ' 



^„(0)z,g;_i(Z;)-(^„-l(0)4;(z,) 
ZyK„4;_j(Z;) 



Proof. These six formulae follow from the Corollaries 13.31 and 13.41 and the bilinear 
product expressions - in particular (|5.6b is a consequence of Il3.87b and Il4.16l4.18l l, 
dSTl from (|3J8j and Il4.16l4.18b . dSlHIl from (|3J9j and J4.17l4.23b, (jSJll is a conse- 
quence of l[3S5t and (I4.17l4.20b . l ISJOll from l l3l96ll and j4.16l4.22b and (|5ll] | from 
(13:971 and (I4.17l4.20b . □ 
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Proposition 5.2. Subject to the conditions in Proposition \5. 1 \ the inverses of the Schlesinger 
matrices (15. 3P and (I5.4I > are given by 

(5.12) 

(0) 

(/>«+! (Z;) 

Kfi Kn+1 ^ 

(pn{Zj) (pn+l{Zj) , 



(5.13) 

( 

(Rr(z;p;))"'=Rr(z;p;-i) = -^^ ■ 



Proof. To establish the first equality of | |5.12| | starting from | |5.3l l and | |5.4| | we require 
(z; Py + 1) and this is found from | |5.4| | under the mapping pj i— > P; + 1. This we 
compute to be the second equality of | |5.12| | and first then follows easily. For the 
first equality of ( 15.131 1 we now require Ri^(z; P; - 1) which is found from ( 15.31 1 under 
the mapping pj i-^ Pj ^ ^- Our computation of this yields the second equality of 
(|5.131 and the first follows in a straight forward way. □ 



(E}:{z-pi))' = -pj;{z;pj + l) = 



(pniZj) 



Kn+1 



(pn+l{Zj) ' 
K„+i(p„{Zj) 



Kn£,n{Zj) 

MO) 

K Z 



C„{Zj) 



-Zi 



From the explicit forms of the Schlesinger transformations it is not obvious that 
they are commutative, however this must be the case and the following proposition 
verifies this fact. 

Proposition 5.3. Any two distinct Schlesinger transformations pj i-> pj+e,pk i-^ pk+£' , 
k + i and e = ±1, e' = +1 commute 

(5.14) Ri'ipj, p, + e')Rf ipj, Pk) = Rf (P; + e, pt)Ri:{pj, p,c). 

Proof. We take the upward shift case {e - e' - +1) first. To check ( |5.14| | we have to 

form Rdpj, Pjc + 1) and P^ipj + 1, Pk) using the explicit formula I I5.3I I which entails 
replacing all coefficients and polynomials appearing therein under the mappings 
Pk p;c + 1 and pj i-> py + 1 respectively. We observe that 

/cicrx I i+,!<+\^ (pn{Zk)(l}„+liZj) - (p„{Zj)(j)„+i{Zk) (k+,j+\^ 

(5.15) [k'„ I = K„K„+2- — ;r^-- — r~\ — I = r« 

^ (p„+l{Zk)(p„+2{Zj) - (pn+-l{Zj)(pn+2{Zk) ^ ' 

All components of the difference between the left-hand and right-hand sides of 
(|5.14t vanish trivially except the (1, 1) component. To show this also vanishes we 
have to use the three term recurrence relation 
(5.16) 

K„+l[(/)„+l(0)(/)„+2(z) - (/)„+2(0)(/)„+i(z)] = z[k„+2(/)„+i(0)(/)„+i(z) - K„(/)„+2(0)(/)„(z)]. 

The other cases of mixed up and down shifted or down shifted parameters can be 
proved using the above equality and employing the inverse formulae | |5.12| | and 
(I5l3t . □ 
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Corollary 5.2. The up-shifted bi-orthogonal polynomials and associated functions have 
evaluations at the singular points given by the formulae 



(5.17) <(z,) = 



(pn+l{Zj)(pn{Zk) - (pnjz j)(pn+l{Zk) , . 

Zj-Zk , K^* J 



(5.18) <!>Tizk) = 



Kn+\(pn{Zj) 

I I / \ . ^n-\-l.\^ / mt-r / \ I j 

(p*„+liZj)(p*„{Zk) - (p'„{Zj)(Pl^^iZk) 



(p„+t{0)(p„{zj) 

K'PniZj) ' 



TniZj)), 



k = j 



(5.19) 



^n{Zk) = 



Kn+l(pn{Zj) ['/'"(^i)^«+l(2't) - (pn+l{ZjKn{Zk)] , ^ ^ 7 

<(/)„+! (0) 



(5.20) S;'^'^(z,) = 



[ )C„K„+i(/)„(2^-) /' 

[(^);(zy)5;,i(z,) - (/);,i(zy).s;(z,)] , ^ ; 



I</'«+l(0)c,')„(Zy) 
^^2z«+i, 



I K„(p„{Zj) i 

provided (pn(zj) + for all n e Z>o, where 



k = j 



(5.21) Tr(z,)=^£^ 

k+i k 



Pk (pniZjWniZk) - (p'„{Zj)(p„{Zk) 



Zi - Zk 



X [</>n(Z/)^;(Zfc) + (pl{Zj)^„{Zk)\ . 

The corresponding results for the down-shifted system is 



(5.22) 0';izk) = 



(5.23) OTizk) = 



Kn-iCn-liZj) 
K„ 



[4>n{Zk)E;„{Zj) + C^l{Zk)UZi)\, k + j 



2z", 



k=j 



.K„-ir„_i(Z;) i' 
K „Zj^i^iz,) bl^Ui^m^k) - Zk^U^j)cp;_,{Zk)] , k^j 

k;^»(0) 



K„K„-iC„-i{Zj) i' 
i KZ UzACniZk) - CnizdUiZk) , 



k = j 
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provided £,'„(zj) + Ofor all n e Z>o, where 



(5.26) T'-{zj) = ^y 



pk ^n{Zj)£,n{Zk) - E;„{Zj)E,n{Zk) 



2z« 



■Zk 



X [E,'„{zj)(\)„{zk) + 5„(z^)(/);,(zO] • 



Compatibility constraints need to be satisfied by the over-determined linear 
bi-orthogonal system and we have to consider three such constraints. The first is 
the compatibility of Schlesinger transformations \5.1\ with the recurrence relation 
I I2.20I I and this implies the following proposition. 

Proposition 5.4. The recurrence matrix and the Schlesinger transformation matrices of 
the bi-orthogonal system satisfy the compatibility condition 



(5.27) 



K+i'y^' Pi)^n{z; pj) = K„{z; pj + l)R'n{z; py). 



Proof. Taking the + case first we construct K„{z;pj + 1) using the specialisation 
a Zj of I l3.88l3.89|l . Forming the difference between the left-hand and right-hand 
sides of (|5.27b we find the (1,2) and (2,2) elements vanish trivially. Then employ- 
ing the forward recurrence for (/)*^j(zy) we find the (2, 1) element also vanishes. 
Utilising the forward recurrences for (pn+i{zj) and cpn+iizj) and finally the previous 
recurrence we see that the (1, 1) element is also identically zero. In the — case we 
compute Kn{z;pj - 1) using the specialisation —> Zj in I|3.96l3.97b . Forming the 
difference again we find (1, 1) element vanishes when the backward recurrence 
I I2.2OI 1 is applied to (5*_j(Zy) and £,*„{zj). The (2, 1) element vanished after application 
of the forward recurrence for 4* ^ ^ (zy) and the backward recurrence for _^ (zy) . The 
(1,2) element vanishes after the backward recurrence for £,n{zj) is utilised. Finally 
the (2, 2) element was shown to be zero upon the application of (|2.8| l, the forward 
recurrence of <^*^^(zy) and the backward recurrence for <^*_^(zy). □ 

The second constraint is the compatibility of the Schlesinger transformations 
with the spectral derivative | |2.27|I . Again we need a preliminary result giving 

explicit formulae for the transformed residue spectral matrices A^^^^^ = A„^k{pj ±1), 
which is contained in the following lemma. 

Lemma 5.1. The transformed residue of the spectral matrix under the Schlesinger trans- 
formation pj pj + 1 is given by 

(5.28) = Rit{zk)A,,k (RiHzk))'' , 

for ki= j or explicitly by 

(5.29) 



A^' = 1 

K„K„+i(p„{Zj)(p„+i{Zj) 



( A. ( \^ (pn+\{OW„{Zj) (p„+l{0)(p„{Zj) 
Kn+l(pn[Zj) H 



-Zk 



Zj - Zk 
K,i+l(p'„{Zj) 
Zj - Zk 



■Zk 



K„+lCp„{Zj) 

Zj-Zk 

(pn+l{0)(p„{Zj) 



A„i 



ZjK„+i(p„(Zj) 

ZkK„+i<pl{Zj) (Zj - Zk)K„+i(p„{Zj) + cl)„+i(0)(p;,(Zj)j ' 



42 



N.S. WITTE 



and 



(5.30) A 



K„-iK„Zy<5;_j(Z;)<5;,(Z/) 



Z; - Zt- K„<;„^Zy; + z 



{Zj - Zk)K„£,;,{zj) + Zjci)„(0)£,„{zj) -ZjK„£,„(zj) 

-Zk(pn{0)i'„{Zj) ZjK„^'„(Zj) 



In the special case k = j we have 
(5.31) A 



_ 1 I -(/)„+! (0)(/);(z,) (0)c/)„(zy) 



"■' K„(/)„+i(zy) \ Z;K„+i(/);(z^) -ZjK„+i(j)„{Zj)j 

X frn I 1U I ^ V ^"-fc /z,K„+i(/)„(z,-) (|)„+i(0)(/)„(Zy)\l 

X ^^p^ + + K„,1<|)„(Z,) - Z, iz,K„.l<|';(Z;) (|)„.1(0)C/);(Z,); /• 



and 



(5-32) A^>-i °)(P;-1) 



2 

J / K„.£;,(Z^) K„<S„(Zy) \y /-^„(0).5„(zy) K„^„{Zj) \ 

Kn-lKni*„_^[Zj)^'n[Zj) ^^„(0)«;(Zy) (^„(0)4„(z^)j Zj - Zfc ^ <^„(0)5;(z^) -Kn^'„{Zj)j " 

Again we require that (p„{zj) + for all n e Z>o for the upward transformation and 
^*„{zj) + Ofor all n e Z>o/or the downward shift. 

Proof. We offer two types of proof here. The first utihses the transformed bi- 
orthogonal polynomials and associated functions given by | |5.17H5.2T] | or l l5.22H5.26l l 
in the formula l l4.6t . One then has to separate out from each of the elements of 
the transformed A„ matrix the original matrix and perform a matrix factorisation 
with the structure of a similarity transformation. This yields the final formulae 
l|5. 2915. 3115. 301 15.32|l . Alternatively one can start from the general Christoffel- 
Geronimus-Uvarov transformation of A„ given by (|3.104l l or (|3.106l l, insert the 
partial fraction expansion (|4.51 and resolve the resulting products on the right- 
hand side of the formula in this way. In order to eliminate the unwanted terms on 
the right-hand side we need to note the eigenvector identities 



^"'\^;,{z,))-^\<P^,{Zj))' 

^«'^^-5;(z,)j- P\-E,*,[z^j' 

[(^*„{Zj) -(p„{Zj))A„,j = -pj{(^*„{Zj) -(\)n{Zj)), 
(E,l(zi) 5„(Z;))A„,, = 0(r„(Z;) UZj)), 

and the form for A„^a 



Using the results gained up to this point we are now in a position to deduce the 
compatibility of the Schlesinger transformations with the spectral derivative and 
deformation derivatives. We begin with the spectral derivative. 
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Proposition 5.5. The spectral matrix and the Schlesinger transformation matrices of the 
bi-orthogonal system satisfy the compatibility condition 



(5.33) 



^R'niz; pj) + R^f (z; p;)A„(z; pj) = A„(z; pj ± (z; pj). 



Proof. Take the + case first. Using | |5.3| | we compute the difference between the 
left-hand and right-hand sides of l|5.33b and resolve the expression into partial 
fractions in the z variable. The (z - z^)"^ term has a coefficient proportional to 
(5.34) 

l(P„^i{Omzj) -c/)„+i(0)</)„(zM ;W(|)„+i(0)c/)*(Z;) -(/)„+i(0)c/)„(zM 

\-ZjK„+i(p*„(Zj) ZjKn+l(pn{Zj) )^ ",i \-ZjKn+l(pniZj) ZylC„+ic/)„(zy) j' 

Using l|5.3Hl this simplifies to 

-cp„^i(0)(P'„{zj) cp„^i{0)cp„{zj)\ 



(5.35) 



ZjK„+i(p*„{Zj) -ZjK„+iCt)„(Zj) 



{pj + A„,j), 



which, with the formula | |4.6t , is identically zero. The nett coefficient of the (z -Zj) ^ 
terms is proportional to 



(5.36) 



<P',{zi)l<Pn{Zj) o)^«'^ Ki[<p'^{Zj)l<Pn{zi) 



/c^y=0 



-rn+l(p'„iZj)/(p„{Zj) rn+ 



Zj(p'„(Zj)/(p„{Zj) 



j+ (-rn+l(pn(Zj)/(p„{Zj) ' n+1 ' 



k\ Zj(p'„(Zj)/(P„{Zj) 



Taking the summand of the above expression, substituting I l5.29b and simplifying 
the result and then employing l|4.6l l we arrive at 



(5.37) 



pk 



Izl K„+i(pl{Zj) 



- {Z, - Z,)^„,,U-,)UZk) ^ -<^„(z,.)</';(Z;)J 



+ ((/)„(Zy)<5;,(Zfc) + c/);(z^)<5„(z,c))((/)„(Z;)(/);(z,c) - (/);(Zy)(^„(z,c)) 



-Z,K„+i Ol 



When this is inserted back into (|5.36l l the first term is amenable to l|4.12|l and cancels 
with another term. Next ll5.31b is utilised in the simplified version of l|5.36l l and 
another term is annulled, leaving only two summation terms. Again appealing 
to (|4.6l l we find the remaining two terms are equal and opposite in sign. Thus 
(|5.36l l is zero. Finally we are left with the (z — z^)"^, + j terms whose coefficient is 
proportional to 



( 



(5.38) 



s , (\>n+\iOWn{Zj) (/)„+l(0)c/)n(Zy) ^ 
Cn+l(pn{Zj) + z . _ Z/ - Zfc 



-Zk 



K„+l(p;,iZj) 



Zj - Zk 
Kn+l(pniZj) 



Zj - Zf; 

Kn+\(p„{Zj) + 
-Zk 

and utilising l |5.29|l we find this to vanish. 



n,k 



I Zj - Zk 

(pn+l{0)(p'„{Zj) 



Zj - Zk 
Kn+l(p'„{Zj) 
Zj - Zk 



(p„+l{0)(pn{Zj) 
Zj - Zk 
Kn+l(pn{Zj) 
Zj-Zk 
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We now treat the — case. Employing l l5.4t and computing the difference of the 
left-hand and right-hand sides of I l5.33|l we can resolve it into partial fractions in z. 
The (z - Zj)~^ terms have a coefficient proportional to 



) (/)„(0)4„(Z;)I 



The first term is easily seen to be zero using l l4.6t . In the second term we employ 
(|5.32t and note that the matrix multiplications yield a null result. The coefficient 
of the (z - Zk)~^ terms is proportional to 



^^■^^^ [zkUmni^j) Z;K„_i4;_i(z,)-Z,K„5;(z,)j^«''^ 

-A'~i ^'^"^"^^i^ ZjKn£,„{Zj) \ 

n,k \zk$„{0)£,'„{Zj) ZyK„-i^;_j(Zy) - ZkKn^'„{Zj)) ' 



Substituting for ^ from | |5.30|I and carrying out the matrix multiplication we find 
this expression is identically zero. Lastly we examine the constant term of the 
partial fraction expansion, which simplifies to 



This is easily shown to be zero using | |4.8|I . This concludes the proof. 



The final constraint is the compatibility of the Schlesinger transformations with 
the deformation derivative | |4.26|I which is established in the following proposition. 

Proposition 5.6. The deformation matrix and the Schlesinger transformation matrices of 
the bi-orthogonal system satisfy the compatibility condition 



(5.42) ^R^f (z; pj) + <(z; pj)B„{z; pj) = B„{z; pj ± l)Rf (z; pj). 

Proof. Our strategy is to show that the difference between the left-hand and right- 
hand sides of | |5.42| | is zero. We begin with the + case first. Our initial task is 

to compute j^Rl^{z) starting with the expressions for R^,^(z) given by | |5.3l l, i5.6i 
and employing the derivatives | |4.30| |, | |4.31| l and | |4.33| |. This yields the following 
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complicated expression (omitting an overall factor of K^lxn) 



(5.43) 



'1'n+lPn,j 



{Z-Zjf\ ZjP„,j 

1 1 



-2, 



— j[z — Zj — rn+lPn,j] 



2 - 2y 2; + r„+iP„,y \ ^ \[-52 + Zj + rn+lPn,j\Pn.j -\zj - r„+iP„_j 



L 

k=l 



pk 



Zk 
Zk 



(pn{zk)£.'„{zk) / i r„+i [2 + Zj + rn+iP„j] \ rl 



22" 



-2[2y + \r„+iP„^j\ 



^ n+1 



(p„{zk)£,„{zk) 



2z" 
^^k 



-Zk 



42P2 . 
2 



■E 



Pk 



\[z + 2j + T„+iP„^j]P„j Zj + jr„+iP„j 

Z; - Zk [(jhAZk)Pn,j - <\>'n{ZkW,n{Zk)Pn,] + ^^(Zfc)] 
7" - Zt 



22" 



1,. 



r„+i[z + Zj + rn+iPn,j] \r^, 



11+1 



z[Zy + 2''n+l-fnj] 

From the definition | |4.27|| we can recast B„(2; py) as 



;c=i 



*(zO.£;,(z;t) / 0\ (p„(zk)^'„{zk) l-l; 



22" 



-1 



22" 







E 

ic=l 



A„ 



Zk 



Z-Zk 



1+ 



In regard to the formula for B„{z; Py + 1) we require deformation derivatives of K„ 

and rl^ . We have found the former in the above calculation whereas the latter can 
be computed starting with l|5.81 along with l|4.33t . In summary we find B„(2; P; + 1) 
to be 



(5.45) 



2; + r„ 







n+lPijj \ Pn,j 2 



EZk 
^'Tk 



k=\ 



(p'„{ZkW„{Zk) 



1 



22^ Zj + rn+lP„,j 

^n{ZkW„(Zk)l 1 



) '"n+1 



22" 
^^k 







- E 

X 



(/)n(Zit).^n(2t) 2;tPK,; 

22^ Zj + r„+iP„,y \-P„,y \ 
Zj-Zk l(p„{Zk)Pn,j - (p'„{Zk)][in{Zk)Pn,j + ^'„{Zk)] 



2z" 

^^k 



2jn+\ 
Zj + rn+lPn,j \ -Zj 



/ I n,i 



k=l 



i+ Zk 
^Z-Zk' 
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Having all the terms in l|5.42l l we now perform a partial fraction decomposition 
with respect to z and examine the resulting coefficients in this decomposition. 
Collecting the l/(z — Zk) terms with k i= j we find its nett coefficient is proportional 
to (|5.38l l in the preceding proposition and thus vanishes. Considering the l/(z - zy)^ 
terms we find the resulting coefficient is proportional to (|5.341 and also vanishes 
according to the preceding proposition. Assembling all the constant terms and 
carrying out the matrix multiplications we find the single terms and summands 
cancel internally without recourse to identities. This leaves the l/(z - Zj) terms 
which we will break down into further sub-divisions. Of these the terms with the 
factor Zj/[Zj + r„+iP„^j] (excluding those proportional to the monodromy exponent 

n - po) are easily seen to cancel. In the summations over A„,fc and A^^""^ we make 
the substitution Zk = Zj - {zj - Zk). We then collect the terms with the factor of 
Zj resulting from this and all the remaining terms with this factor. Using l|5.36t 
these terms are simplified to (they do not vanish in totality because the sum over 
k excludes k = Q) 



(5.46) - - 

L 



1 

—TfiZjPfij I 

The A: = term is extracted from the summations containing Zj - Zk and simplifies 
to 

(5-47) |(«-po)(l-r„P„,;)P„.;|':;J oj. 

Adding l |5.46|l and (|5.471 we get 

(5.48) Zj{n - po)rn(^2'' _p^J. 

Now consider the k terms in the summations containing zy - Zk- Using the 
working leading to (|5.471 and the formulae l|5.29b and 1 14.61 we find the summands 
reduce to 

(5.49) {±j-tk){A„,kh,2(-P^'' P^]- 

\ n,j "'I I 

Combining the contributions of (|5.48l l and (|5.491 we have 

(5.50) Yj^^'~^'^^^"''^A-P^'' pI)- 

There is one group of terms in the l/(z - Zj) coefficient to be simplified, namely 
the summation with the factor Zk/Zk- After considerable internal cancellation this 
group simplifies to 

(5.51) YjZkiAn,kK2i^l2' p^J. 
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In summary the l/(z - zj) coefficient reduces to 

(5-52) Zj ^{An,k)l,2 (_p2'''' J oc -Zj{An,cc)l,2 = 0, 

where we have used (|4.8l l in the last step. 

Next we treat the - case. We recall the definition of Q„j := £.„{zj)/^*^{zj). We 

compute ■ffRi't (2; Pj) first beginning with l l5.4t and ||5.9|I and utilising the derivatives 
(|4.321 and l|4.341 . Up to an overall factor of )c„/k~ we find this term is 



(5.53) 



ZL 

2z, 



ZjQn,j 

-f„Z Z + Zy(l + f„Q„,;) 



Jc=l 



(j);,(zk)£.;,{zk) 



2z" 



Z[-1 + 5f„Q„,y] "IQhjIz + Zjil + r„Qn,j)]j 



-E 



Pk-. 



(Pn{Zk)C„{Zk) I -z^Q„,y , 
2z^ \zf„ 

Zy - Zfc [(/)„(zO + Q«j</';(Zfc)][Q«j5*,(Zic) - ^n{Zk)] 



k^j=0 



■Zk 



2z" 



2Z j^n 

1 -2 

-r 7 

2 n 



1 + J'nQn,; 

jf„[z + Zj{l+r„Q„_j)]j 



However in contrast to the formula for B„(z; pj) (|5.44l l the one for B„(z; p; - 1) is 
more complicated as we have to compute the deformation derivatives of the shifted 
quantities k~, f^. Starting with the expressions | |5.9| | and | |5.11| | we can employ the 
derivatives computed earlier in addition to | |4.30| |. In summary we find B„(z; py — 1) 
to be given by 



(5.54) ^ i 



r„ - 



L 

k=l 



pk 



Zk 
Zk 



[(p'„{Zk) - T„(p„{Zk)][£.''„{Zk) + fnlniZk)] 



Zk 



2z" 



(\>l{zk)E,l{zk) I -lQ„^j 



2z" 
^^k 



-[1 + rnQn,j] iQnJ ^ 



Z;[l + f„Q„,y] 
(p„{Zk)i;,{Zk) I -I 



0\ 

1 



2z" 







E 

k*j=0 



PK 



■ Zk [(p„{Zk) + Qn,j(p*„iZk)][Q„,j£,;,{Zk) - ^n{Zk)] 



■Zk 



2z" 



1 + r„Q„,j 

/ I n,i 



r„ - 



Zk 



k=l 



■Zk 



We now decompose ||5.42| | into partial fractions with respect to z. Taking the 
l/(z - Zk) terms we find they are proportional to the expression (|5.40l l which was 
shown earlier to vanish. Collecting the z terms together one finds they vanish 
identically in a straight forward manner. This leaves the constant terms. After 
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making the obvious simplifications here one is left with 
(5.55) z,ll \]Z^njA-Y^Al-Jl \ 

\ ' \ ' k=l k=l \ 



"ly . [^'ni^k) + r„^„{zk)][(p'„{zk) - f»(/)»(zt)] 

l + f„Q„,;\l QnJlt^'"' 2Z« 

-1 \ y Zj - Zk [(pnjZk) + Qn,j(p'„{Zk)][Q „,j£,'„{Zk) ~ ^n{Zk)] 



In the second and third terms we substitute Zk - Zj - {zj - Zk) and split them into 
two terms each. Collecting the terms with a factor of Zj we can perform the k 
summation using l|4.8l l. This group of terms vanishes because of l l5.4Hl . In the 
remaining group we employ ||5.30|I and | |5.32|I along with | |4.6|I and carry out the 
matrix multiplications yielding 

V f . ,mzk) + rnUzk)]l \ 

(5.56) -yk{z,-Zk) ^^^(^^^j 

^/ _ " [(p,r{^k) + Qn,j(p'„(Zk)][Qn,j^ni^k) " <^»(Z)c)] /f„ -1 \ 



E. , [(pni^k) + Q„,j(t)'„{zk)M'„{zk) + r„^„{zk)] 1 
Pkizj - Zk) ^ 

le=r\ k \ 



l + f„Q,w^^^^^^ 2z« \rn -Ij- 

The second term of the above expression cancels with the last term of ( |5.55t . In 
what remains of the constant terms we separate them into those with a factor of Zj 
and those with Zk- In the former set of terms we can perform the k summation using 
I|4.12ll4.15l l and find complete cancellation. In the latter set we find the summands 
mutually cancel. Thus all the constant terms vanish identically and the theorem is 
proved. □ 

It is expected that the explicit Schlesinger transformations derived here | |5.3| |, 
(|5l4t , ll5l29l l, llESB/ l|530t , | |532|| should co-incide with those foimd by the Kyoto 
school [32J however we refrain from carrying out this verification as this would 
entail a significant enlargement of the tasks considered here and add substantially 
to the length of the present study. 

We conclude this section and our study with a selection of results on the bi- 
linear difference equations satisfied by the T-functions of our integrable system, i.e. 
the Toeplitz determinants. These are pure difference equations in the parameters of 
the system, n, {p;};=o,...,M or the monodromy exponents {0;};=o,(...,m,'=o/ products of 
T-functions, and are known as Hirota-Miwa equations. So in this sense we return 
to one of the initial motivations of our study, namely to find alternative means of 
characterising averages over the unitary group (|1.2t of the regular semi-classical 
weights | |4.3|| . All of the Hirota-Miwa type equations are very easily derived from 
the preceding analysis and one might consider them to be immediately apparent. 
Let us define the following notation for the Toeplitz determinants 

(5.57) I„[{pj}] :=t(0o 0;,...,0oc). 
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Proposition 5.7. For j = 1,...,M (the j = cases are trivial) the z-functions satisfy the 
Hirota-Miwa equations 

(5.58) t()t(0o + 1, 0OC + 1) 

= t(0; + 1, 0OO + 1)t(0o + 1, 0; - 1) - ZjT{dj - 1, 0OO + 1)t(0o + 1, 0y + 1), 

(5.59) t()t(0o - 1, + 1) 

= t(0; - 1, 0OO + 1)t(0o - 1, 0; + 1) + Z;T(0y + 1, 0oo + 1)t(0o " 1, Oj - 1), 

i(;/!ere t() denotes the T-function with unshifted monodromy exponents. These constitute 
all of the independent equations involving one generic exponent dj. 

Proof. These follow from identifying the bi-orthogonal pol3momials and associ- 
ated functions in terms of modified Toeplitz determinants using I l3.1t -l l3.4|l and 
employing these formulae in the Casoratian equations l|2.24|l - l|2.26|l . We find that 

T T _ fi+fh _ 7 7/+r/- 

\ I 'ni„ - i-n '-n + ^/-'„_i-'„+i' 
\ J 'ni-n ~ '-n -'n ^ ^;-'„_i-'„+i' 

and we note that these last two equations lead to the same bi-linear equation. 
Consequently we deduce | |5.58l l and | |5.59| |. □ 

Proposition 5.8. For any j + k = 0, . . . ,M the z-functions satisfy the Hirota-Miwa 
equations 

(5.60) (Zy - zOt(0o + 1, Sco + l)T(e; - 1, 0,c - 1, Soo + 2) 

= t(0o + 1, - 1, 0OC + 2)t(0; - 1, 0OC + 1) - t(0o + 1, 0; - 1, 0OO + 2)T(0fc - 1, + 1), 

(5.61) t(0o + 1, 0OO + 1)t(0; - 1, Si + 1) 

= -ZtT(0o + 1, 0JC + 1)t(0; - 1, 0OO + 1) + T(0t + 1, 0OO + 1)t(0o + 1, 0; " 1), 

and 

(5.62) (zy - zOT()T(eo + 1, 0; -hdk- 1, 0o„ + 1) 

= ZjTidj - 1, 0CO + 1)t(0o + 1, 0/c - 1) - 2)cT(0O + 1, 0; " 1)t(0;c " 1, 0o<, + 1), 

and 

(5.63) T()T(eo + 1, 0; - 1, 0t + 1, 0OO - 1) 

= ZjTidj - 1, 0OC - 1)t(0o + 1, 0t + 1) + t(0o + 1, dj - l)T(0t + 1, 00. - 1). 

These constitute all of the independent equations involving two generic exponents dj, dt- 

Proof. Equation (|5.60t is derived by combining the expression for the double shift 
Pj pj + l,pk p/t + 1 given in (|5.15t with the single shifts pj^k ^ Pj,k + 1 as in 
i5.6\ and this results in 

(5.64) - ce,] = i„,,ii::i: [^./r - ca^] 
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This can rewritten as a perfect difference equation and solved, yielding 

i„ l„+i - 1„ Jn+i - (Zj - Zk)ln+U„ 

A similar method applies when deriving the other equations and we confine our- 
selves to recording the respective intermediate steps, namely 

j*j+„k- _ ji+jk- _T Ti+^- 

and 
and 

□ 



6. Acknowledgments 

This research has been supported by the Australian Research Council. The 
author wishes to acknowledge the generous advice given by Peter Forrester, the 
interest shown by Mourad Ismail and the detailed attention given to the work by 
an anonymous referee. 

References 

1. M. Adler and P. van Moerbeke, Recursion relations for unitary integrals, combinatorics and the Toeplitz 
lattice, Comm. Math. Phys. 237 (2003), no. 3, 397-440. MR 1 993 333 

2. G. Akemann and G. Vernizzi, Characteristic polynomials of complex random matrix models, Nuclear 
Phys. B 660 (2003), no. 3, 532-556. MR MR1982915 (2004h:82041) 

3. J. Baik and E. M. Rains, Algebraic aspects of increasing subsequences, Duke Math. J. 109 (2001), no. 1, 
1-65. MR 1 844 203 

4. G. Baxter, Polynomials defined by a difference system, J. Math. Anal. Appl. 2 (1961), 223-263. MR 23 
#A3421 

5. M. C. Bergere, Biorthogonal Polynomials for Potentials of two Variables and External Sources at the 
Denominator, hep-th/S484126, 2004. 

6. M. Bertola and M. Gekhtman, Biorthogonal Laurent polynomials, Toplitz determinants, minimal Toda 
orbits and isomonodromic tan functions, Constr Approx. 26 (2007), no. 3, 383^30. MR MR2335689 
(2008h:42043) 

7. M. I. Bueno and F. Marcell^, Darboux transformation and perturbation of linear functionals. Linear 
Algebra Appl. 384 (2004), 215-242. MR MR2055354 (2005b:15027) 

8. E. B. Christoffel, liber die Gaussische Quadratur und eine Verallgemeinerung derselben, J. Reine Angew. 
Math. 55 (1858), 61-82. 

9. J. B. Conrey, D. F. Farmer, and M. J. Zirnbauer, Hozve pairs, supersymmetry, and ratios of random 
characteristic polynomials for the unitary groups Ui^, arXiv: ma th-ph/85 11824, 2005. 

10. , Autocorrelation of ratios of L-f unctions, arXiv: math. NT/8711 .0718,2007. 

11. J. B. Conrey, P. J. Forrester, and N. C. Snaith, Averages of ratios of characteristic polynomials for the 
compact classical groups, hit. Math. Res. Not. (2005), no. 7, 397-431. MR MR2130839 (2005m:11167) 

12. P. J. Forrester, Log Gases and Random Matrices, 

http : //www . ms . unimelb . edu . au/matp j f /matp j £ . html. 

13. , Growth models, random matrices and Painleve transcendents, Nonl. 16 (2003), no. 6, R27-R49. 

14. P. J. Forrester, N. E. Frankel, T. M. Garonl, and N. S. Witte, Painleve transcendent evaluations of finite 
system density matrices for Id impenetrable bosons. Comm. Math. Phys. 238 (2003), no. 1-2, 257-285. 
MR MR1990877 (2005c:82011) 

15. P. J. Forrester and E. M. Rains, Symmetrized models of last passage percolation and non-intersecting lattice 
paths, arXiv:math-ph/S7S5 . 3925, 2007. 

16. P. J. Forrester and N. S. Witte, Application of the z-function theory of Painleve equations to random matrices: 
Pvi, the JUL, CyUE, cJUE and scaled limits, Nagoya Math. J. 174 (2004), 29-114. MR MR2066104 
(2005h:33047) 



51 



17. , Bi-orthogonal Polynomials on the Unit Circle, regular semi-classical Weights and Integrable Sys- 
tems, Construct. Approx. 24 (2006), no. 2, 201-237. 

18. R. Gamier, Sur des equations differentielles du troisieme ordre dont I'integrale generale est uniforme et sur 
une classe d' equations nouvelles d'ordre superieur dont I'integrale generale a ses points critiques fixes, Ann. 
Sci. Ecole Norm. Sup. (3) 29 (1912), 1-126. MRMR1509146 

19. Ya. L. Geronimus, Orthogonal polynomials: Estimates, asymptotic formulas, and series of polynomials 
orthogonal on the unit circle and on an interval. Authorized translation from the Russian, Consultants 
Bureau, New York, 1961. MR 24 #A3469 

20. E. Godoy and F. Marcellan, Rational modifications of Borel measures. Proceedings of the Xllth 
Portuguese-Spanish Conference on Mathematics, Vol. II (Portuguese) (Braga, 1987) (Braga), Univ. 
Minho, 1987, pp. 326-330. MR MR1139178 

21. , An analog of the Christoffel formula for polynomial modification of a measure on the unit circle, 

BoU. Un. Mat. Ital. A (7) 5 (1991), no. 1, 1-12. MR MR1101005 (92c:42022) 

22. , Orthogonal polynomials and rational modifications of measures, Canad. J. Math. 45 (1993), no. 5, 

930-943. MRMR1239908 (95a:42031) 

23. J. Hamad and A. Yu. Orlov, Integrals of rational symmetric functions, Two-Matrix Models and Bi- 
orthogonal Polynomials, math-ph/9693949, 2006. 

24. , Fermionic approach to the evaluation of integrals of rational symmetric functions, 

math-ph/6794 . 1159vl, 2007. 

25. E. Hendriksen and H. van Rossum, Orthogonal Laurent polynomials, Nederl. Akad. Wetensch. Indag. 
Math. 48 (1986), no. 1, 17-36. MR MR834317 (87i:30008) 

26. M. E. H. Ismail and R. W. Ruedemann, Relation between polynomials orthogonal on the unit circle with 
respect to different weights, J. Approx. Theory 71 (1992), no. 1, 39-60. MR 941:42028 

27. K. Iwasaki, H. Kimiua, S. Shimomura, and M. Yoshida, Prom Gauss to Painleve, Aspects of Math- 
ematics, E16, Friedr. Vieweg & Sohn, Braimschweig, 1991, A modern theory of special fimctions. 
MRMR1118604 (92j:33001) 

28. M. Jimbo and T. Miwa, Deformation of linear ordinary differential equations. I, Proc. Japan Acad. Ser. A 
Math. Sci. 56 (1980), no. 4, 143-148. MRMR575994 (84k:14009) 

29. , Deformation of linear ordinary differential equations. II, Proc. Japan Acad. Ser. A Math. Sci. 56 

(1980), no. 4, 149-153. MR MR575994 (84k:14009) 

30. , Deformation of linear ordinary differential equations. Ill, Proc. Japan Acad. Ser. A Math. Sci. 56 

(1980), no. 6, 269-274. MR MR581469 (85i:58005) 

31. , Studies on holonomic quantum fields. XVII, Proc. Japan Acad. Ser. A Math. Sci. 56 (1980), no. 9, 

405-410. MR 85h:82016 

32. , Monodromy preserving deformation of linear ordinary differential equations with rational coeffi- 
cients. 11, Phys. D 2 (1981), no. 3, 407-448. MR 83k:34010b 

33. M. Jimbo, T. Miwa, and K. Ueno, Monodromy preserving deformation of linear ordinary differential 
equations with rational coefficients. 1. General theory and t-furwtion, Phys. D 2 (1981), no. 2, 306-352. 

MR 83k:34010a 

34. W. B. Jones and W. J. Thron, Survey of continued fraction methods of solving moment problems and 
related topics, Analytic theory of continued fractions (Loen, 1981), Lecture Notes in Math., vol. 932, 
Springer, Berlin, 1982, pp. 4-37. MR MR690451 (84k:44015) 

35. W. B. Jones, W. J. Thron, and O. Njastad, Orthogonal Laurent polynomials and the strong Hamburger 
moment problem, J. Math. Anal. Appl. 98 (1984), no. 2, 528-554. MR MR730525 (85h:44015) 

36. H. Kimura, Symmetries of the Gamier system and of the associated polynomial Hamiltonian system, Proc. 
Japan Acad. Ser. A Math. Sci. 66 (1990), no. 7, 176-178. MR 91j:58060 

37. A. P. Magnus, Painleve-type differential equations for the recurrence coefficients of semi-classical orthogonal 
polynomials. Proceedings of the Fourth International Symposium on Orthogonal Polynomials and 
their Applications (Evian-Les-Bains, 1992), vol. 57, 1995, pp. 215-237. MR 961:42027 

38. M. Mazzocco, The geometry of the classical solutions of the Gamier systems. Int. Math. Res. Not. (2002), 
no. 12, 613-646. MRMR1892536 (2002m:34137) 

39. B. McCoy and T. T. Wu, The Two-Dimensional Ising Model, Harvard University Press, Harvard, 1973. 

40. K. Okamoto, Studies on the Painleve equations. I. Sixth Painleve equation Pyi, Arm. Mat. Pura Appl. (4) 
146 (1987), 337-381. MR 88m:58062 

41. P. I. Pastro, Orthogonal polynomials and some q-beta integrals of Ramanujan, J. Math. Anal. Appl. 112 
(1985), no. 2, 517-540. MR MR813618 (87c:33015) 

42. R. W. Ruedemann, Some results on biorthogonal polynomials. Internal. J. Math. Math. Sci. 17 (1994), 
no. 4, 625-636. MR MR1298784 (95j:33026) 



52 



N.S.WnTE 



43. L. Schlesinger, Uber eine Masse von Dijferentialsystemen beliebiger Ordnung mitfesten kritischen Punkten, 
J. Reine Angew. Math. 141 (1912), 96-1145. 

44. B. Simon, Orthogonal polynomials on the unit circle. Part 1, American Mathematical Society Collo- 
quium Publications, vol. 54, American Mathematical Society, Providence, RI, 2005, Classical theory. 
MRMR2105088 (2006a:42002a) 

45. , Orthogonal polynomials on the unit circle. Part 2, American Mathematical Society Collo- 
quium Publications, vol. 54, American Mathematical Society, Providence, RI, 2005, Spectral theory. 
MR MR2105089 (2006a:42002b) 

46. V. Spiridonov, L. Vinet, and A. Zhedanov, Spectral transformations, self-similar reductions and orthogonal 
polynomials, J. Phys. A 30 (1997), no. 21, 7621-7637. MR MR1603364 

47. T. Suzuki, Affine Weyl group symmetry of the Gamier system, Fimkcial. Ekvac. 48 (2005), no. 2, 203-230. 
MRMR2177118 (2006h:33014) 

48. G. Szego, Orthogonal polynomials, third ed.. Colloquium Publications 23, American Mathematical 
Society, Providence, Rhode Island, 1967. 

49. T. Tsuda, Birational symmetries, Hirota bilinear forms and special solutions of the Gamier systems in 
2-variables, J. Math. Sci. Univ. Tokyo 10 (2003), no. 2, 355-371. MR 1 987 136 

50. V. B. Uvarov, Relation between polynomials orthogonal with different weights, Dokl. Akad. Nauk SSSR 
126 (1959), 33-36. MR 26 #6679 

51. , The connection between systems of polynomials that are orthogonal with respect to different distri- 
bution functions, USSR Comput. Math, and Math. Phys. 9 (1969), 25-36. MR 41 #7369 

52. N. S. Witte, Bi-orthogonal systems on the unit circle, regular semi-classical weights and the discrete Gamier 
equations., 2007. 

53. A. Zhedanov, Rational spectral transformations and orthogonal polynomials, J. Comput. Appl. Math. 85 
(1997), no. 1, 67-86. MR MR1482157 (98h:42026) 

54. , The "classical" Laurent biorthogonal polynomials, J. Comput. Appl. Math. 98 (1998), no. 1, 

121-147. MRMR1656982 (99k:33065) 

Department of Mathematics and Statistics, University of Melbourne, Victoria 3010, Australia 
E-mail address: n . witte@ms . unimelb . edu . au 



